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1. Introduction: the Heisenberg-Euler Effective Lagrangian 

1.1. The Spinor and Scalar QED one loop results 

In classical field theory the Lagrangian encapsulates the relevant classical 
equations of motion and the symmetries of the system. In quantum field 
theory the effective Lagrangian encodes quantum corrections to the classical 
Lagrangian, corrections that are induced by quantum effects such as vac- 
uum polarization. This can be used as a semi-phenomenological device, as 
in effective field theory, or as a fundamental approach in which one uses 
an external classical field as a direct probe of the vacuum structure of the 
quantum theory. The seminal work of Heisenberg and Euler [1], and Weis- 
skopf [2] produced the paradigm for the entire field of effective Lagrangians 
by computing the nonperturbative, renormalized, one-loop effective action 
for quantum electrodynamics (QED) in a classical electromagnetic back- 
ground of constant field strength. This special soluble case of a constant 
field strength leads immediately to several important insights and applica- 
tions. 

In spinor QED, the one-loop effective action for electrons in the presence 
of a background electromagnetic field is 

S'(i) = -ilndet{il/)-m) = -- lndet(^ + m^) (1.1) 

where the Dirac operator is ^= 7*^ (d^, + ieA^), Ai, is a fixed classical gauge 
potential with field strength tensor F^j^ = d^Ay — dyA^^ and m is the electron 
mass. This one-loop effective action has a natural perturbative expansion 
in powers of the external photon field A^, as illustrated diagrammatically in 
Figure 1. By Furry's theorem (charge conjugation symmetry of QED), the 
expansion is in terms of even numbers of external photon lines. Heisenberg 






Figure 1. The diagrammatic perturbative expansion of the one loop eflfective action (1.1). 

and Euler, and Weisskopf, showed that in the low energy limit for the ex- 
ternal photon lines, in which case the background field strength F^^ could 
be taken to be constant, it is possible to compute a relatively simple closed- 
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form expression for the effective action, which generates all the perturbative 
diagrams in Figure 1. Heisenberg and Euler [1] expressed their final answer 
for spinor QED in several equivalent ways: 
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Here £c is the critical field strength 
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and a and h are related to the Lorentz invariants'^ characterizing the back- 
ground electromagnetic field strength [1]: 



-h^ = E^ -B' 
ab = E-B = 



FnuF^" = -2T , 






-Q ■ 



Thus 



^j:2 +g2_j: 



^/J^+G^+T . 



:i.4) 
:i.5) 

:i.6) 



If ^ 7^ 0, it is possible to transform to a Lorentz frame in which the elec- 
tric and magnetic fields are parallel or antiparallel, depending on the sign 



''There has been a notational reversal [3,8,62] of a <-> 6 since the original Heisenberg-Euler 
paper [1]. I stick here with Heisenberg's original notation since in a frame in which B and E are 
parallel, we associate b '^ B and a *^ E, which seems more natural. Also, modern formulations 

2 2 

have adopted Schwinger's choice [12] of units in which the fine structure constant a = j- * f~i 

in which case the prefactor in (1.2) is -r > o"?- 
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oi Q. A suitable sign choice is implicit in (1.6). The a and b parameters 
are significant because ±6 and itia are the eigenvalues of the 4x4 antisym- 
metric constant matrix F^jy. This explains why this constant field strength 
case is exactly soluble; a constant F^i, can be represented by a gauge field 
Af^ = —^Ff^^x'^, which is linear in x. Thus, in an appropriate basis, the 
Dirac operator factorizes into two independent Landau level problems, of 
"cyclotron" frequencies b and ia respectively. Hence the traces in (1.1) can 
be done in closed form in terms of trigonometric functions, leading directly 
to (1.2). For details see [1-8]. 

Weisskopf [2] computed the analogous quantity for scalar QED 



5£L = |lndet(D2+^2 



:i.7) 



which involves the Klein-Gordon operator rather than the Dirac operator. 
Weisskopf obtained (I have converted from his different units and different 
definition of the critical field, which introduces various factors of 27r): 
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Note that the prefactor becomes t^-^ in Schwinger's units [12]. 



1.2. Physical Applications 

The Heisenberg and Euler result (1.2) leads immediately to a number of 
important physical insights and applications. 
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1.2.1. Nonlinear QED Processes : Light- Light scattering 

The Euler-Heisenberg effective Lagrangian (1.2) is nonlinear in the elec- 
tromagnetic fields, the quartic and higher terms representing new nonlinear 
interactions, which do not occur in the tree level Maxwell action. The first of 
these new interactions is light-light scattering, represented diagrammatically 
by the second Feynman diagram in the expansion in Figure 1. Expanding 
the Euler-Heisenberg answer to quartic order we find 



e 



4 



5(1) = - / d^x 



{E^ - B^Y + 7{E-BY +... (1.9) 



which gives the low energy limit (since the field strength was constant) of the 
amplitude for light-light scattering. As first discussed by Euler and Kockel 
[9], these nonlinearities can be viewed as dielectric effects, with the quantum 
vacuum behaving as a polarizable medium. In Weisskopf's words [2]: 

"WLen passing through electromagnetic fields, light will behave as if 
the vacuum, under the action of the fields, were to acquire a dielectric 
constant different from unity." 

The full light-light scattering process in QED was not solved until 1951 by 
Karplus and Neuman [10]. 

1.2.2. Pair-production from vacuum in an electric field. 

The presence of a background electric field accelerates and splits virtual 
vacuum dipole pairs, leading to e^e~ particle production: see Figure 2. 
This instability of the vacuum was realized already by Heisenberg and Euler 
[1], motivated in part by earlier work of Sauter on the Klein paradox [11]. 
This pair production process was later formalized in the language of QED 
by Schwinger [12, 13]. Heisenberg and Euler [1] deduced the leading pair 
production rate in a weak electric field to be 



47r3 



exp 






(1.10) 



This rate is deduced from the imaginary part of the effective Lagrangian 
(1.2) when the background is purely electric 

r = 2/m£. (1.11) 

In modern language, this imaginary part gives the rate of vacuum non- 
persistence due to pair production [12,13]. The rate is extremely small for 
typical electric field strengths, becoming more appreciable when the E field 
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Figure 2. A static electric field can tear apart a virtual e+e^ pair from the vacuum, producing 
an asymptotic electron and positron, as shown on the left. On the other hand, a static magnetic 
field does not break this virtual dipole apart, as shown on the right for a magnetic field directed 
out of the page. 

approaches a critical value Ec ~ ™^ ~ 10^^ Vcm~ , where the work done 
accelerating a virtual pair apart by a Compton wavelength is of the order 
of the rest mass energy for the pair. Such electric field strengths are well 
beyond current technological capabilities, even in the most intense lasers. 
For an excellent recent review of the search for this remarkable phenomenon 
of vacuum pair production, see [14]. Even though the condition of a con- 
stant electric field is rather unrealistic, Heisenberg and Euler's result (1.10) 
provides the starting point for more detailed analyses which incorporate 
time-dependent electric fields, as is discussed below in Section 2. 



1.2.3. Charge renormalization, f3-functions and the strong-field limit. 

Another remarkable thing about Heisenberg and Euler's result (1.2) is that 
they correctly anticipated charge renormalization. The first term (on each 
line) on the the RHS of (1.2) is the bare result, the second term is the 
subtraction of a field-free infinite term, and the third term is the subtraction 
of a logarithmically divergent term which has the same form as the classical 
Maxwell Lagrangian. This last subtraction corresponds precisely to what 
we now call charge renormalization, as was later formalized by Schwinger 
[12, 13]. Indeed, the study of such logarithmically divergent terms was a 
major focus of the early quantum field theory work of both Heisenberg and 
Weisskopf. Weisskopf [2] noted the characteristic strong-field limit behavior 
of the Heisenberg-Euler result (1.2), for example for spinor QED in a strong 
magnetic background: 

^^P"'^"-' - '' log (^) , B^oo. (1.12) 



-^Maxwell 127r2 V 



m 



In modern language, the coefficient of the logarithmic dependence of this 
ratio is known as the one-loop QED /^-function, and Weisskopf anticipated 
the importance of such logarithmic behavior. In later work [15] he showed 
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that for n > 2 foop order in perturbation theory, the divergence was at most 
log"~ , a fact that was important for the work of Gell-Mann and Low [16] and 
the devefopment of the renormahzation group. This connection between the 
strong field Hmit of effective Lagrangians and /3-functions wih be discussed 
in more detail below in Section 5. 

1.2.4. Low-energy effective field theory 

The Heisenberg-Euler result (1.2) is the paradigm of what is now called "low 
energy effective field theory" [17, 18]. In this approach one describes the 
physics of some light degrees of freedom (here the photon field) at energies 
much lower than some energy scale above which one has integrated out the 
heavy degrees of freedom (here the electron field) . Generically, the effective 
Lagrangian is expanded in terms of gauge and Lorentz invariant operators 
QW for the light fields, respecting the relevant remaining symmetries: 

^cS = m^Y,an—;^ . (1.13) 

n 

By power counting, the operators O^"' are balanced by appropriate pow- 
ers of the heavy mass scale m. We can see this structure directly in the 
Heisenberg-Euler result (1.2), for example in the first nontrivial term, the 
light-light scattering term (1.9). The full expansion is given below in (1.15) 
and (1.17). The light photon field is described in terms of gauge and Lorentz 
invariant operators constructed from the field strength F^jy , and having mass 
dimension n. At mass dimension 8, we can have {F^yF^^'^)'^ or {F^^F'^'^)'^, as 
in (1.9). At mass dimension 10, we could have terms involving derivatives 
of the field strength, such as {d^F^pd'^F'^P){Fai3F"^), which do not show 
up in the Heisenberg-Euler result, but which can be found in a derivative 
expansion about the constant field result, as discussed below in section 2. 
The effective field theory interpretation of the Heisenberg-Euler Lagrangian 
also makes connection with the operator-product-expansion (OPE), where 
the polarization tensor is expanded as 

n^. = (q^qu - q^g^u) J^ Cn{Q^) (0(")) (1.14) 

n 

as is discussed in detail in [19-21]. 

1.3. Weak-field expansions of Heisenberg-Euler 

In this section I present some results for various weak-field expansions of the 
Heisenberg-Euler spinor QED effective Lagrangian (1.2), as well as for the 
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corresponding scalar QED effective Lagrangian (1.8) derived by Weisskopf. 

1.3.1. Spinor QED case 

The weak field expansion of (1.2), expressed in terms of the Lorentz invari- 
ants a and b defined in (1.5), is: 

^(1) "i^ ST^/r, qm Y^ ^2k^2n-2k /^2e6\ " f 2iea\ 

^spinor ~ -^ Z>( ^ " ^^ Z> (2k)\{2n - 2k)\ y^J \~^) ■ 

(1.15) 

Here the B2n are the Bernoulli numbers [22-25], which arise from the Taylor 
expansions of the trigonometric functions appearing in (1.2): 

°^ o2fcH , °^ (_^\k'^2k) 



k=0 ^ '' k=0 ^ '' 

At any given order in m"^", the combinations of a and h can be rewritten in 
terms of the more familiar Lorentz invariants {E"^ — B'^) and E ■ B m (1.5). 
For example, the first two orders yield: 

(1) ^ \360 "I" 72 "I" 360^ ^ \^630 "•" 180 180 630^ 

-^spinor ;^^4^^2 + ;;^8^^2 + ' ' ' 

_ e* [(a2 - 62)2 + 7(^j,)2] ^6(^2 _ ^2) [2(a2 - fe2)2 + 13(^5)2] 

360m%2 + 1260m%2 "^ ' ' ' ' 

(1.17) 

The first term in (1.17) corresponds to the light-light scattering result (1.9), 
while the second gives the low energy limit of the diagram in Figure 1 with 
six external photon legs, and agrees with Eqn. (43) in [1]. 

Several special cases of the general expansion (1.15) are of interest. 

(1) Purely magnetic background : If the background is purely magnetic, 
of strength B > 0, then the integral representation (1.2) reduces to 

which has the asymptotic expansion 

£« r^-H^y ^11^ ^] (119) 

spinor 87r2^ (2n + 4)(2n + 3)(2n + 2) V"^V 
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This example illustrates clearly that the Heisenberg-Euler perturbative 
weak field expansion (1.15) is a divergent asymptotic series, since the 
Bernoulli numbers B2n alternate in sign and diverge factorially fast in 
magnitude [22,23]: 

^2n = (-ir^2|j^C(2n). (1.20) 

Here C(^) denotes the Riemann zeta function [23] 

oo 

CH = Efc^ (1-21) 

k=l 

which is exponentially close to 1 for large integers n. In this magnetic 
field case, the perturbative series is divergent but alternating, and 
is Borel summable [26]. In fact, the Borel sum of (1.19) is simply 
the original Heisenberg-Euler integral representation (1.18), as can be 
checked using the relation [22,23] 

1 s °° -25^ 

-tl^^ - 7 - 3 = ^ Pvr^(.^ + W) • ('•'') 

fc=i ^ ^ 

(2) Purely electric background : If the background is purely electric, of 
strength E > 0, the integral representation (1.2) reduces formally to 

<s.„^-S^X'°f("-M)-"'""^'- <'-> 

This integral representation has the asymptotic expansion 

,(1) ^^V (-l)"g2.+4 r'^^V"^\l211 

spinor 87r2 ^ (2n + 4)(2n + 3)(2n + 2) V "i2 y -V- 7 



n=0 

If only one oi E or B is nonzero, the only Lorentz invariant is E"^ — B^; 
therefore, changing from the perturbative series (1.19) to the perturba- 
tive series (1.24) simply involves the replacement B^ -^ —E"^. There- 
fore, in the electric case the perturbative series (1.24) is divergent and 
nonalternating, and thus is not Borel summable [26]. This has a dra- 
matic physical consequence. It means that the effective Lagrangian has 
an exponentially small nonperturbative contribution which is imagi- 
nary, even though it is real to all orders in perturbation theory. This is 
a classic signal of instability, and corresponds precisely to the vacuum 
instability of vacuum pair production first noted by Heisenberg and 
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Euler. This imaginary part can also be seen to arise from the poles on 
the real s axis encountered in the integral representation (1.23), which 
shows that this integral representation must be defined more carefully, 
with an approriate ie prescription [3,12]. From this integral represen- 
tation, or from a Borel dispersion relation, it is straightforward to de- 
rive the full nonperturbative imaginary part of the Heisenberg-Euler 
effective Lagrangian to be: 



2 j-,2 °^ 1 
n=l 



m?"Kn 
eE 



(1.25) 



In modern language this result, first derived by Schwinger [12], gives 
the higher order instanton contributions to the leading imaginary part 
(1.10) found by Heisenberg and Euler. 

(3) Perpendicular magnetic and electric background fields : If i? • i? = 0, 
then ^ = 0, and from (1.6) we see that either a = or 6 = 0, depending 
on the sign of J-. If E"^ > B^, then the situation is as for the purely 
electric field, with a vacuum instability. On the other hand, if B^ > 
E"^, then the system behaves like the case of a purely magnetic field. If 
E"^ = B^, then jC], \^^^ vanishes and there are no quantum corrections 
at all. 

(4) Parallel magnetic and electric background fields : As mentioned pre- 
viously, ii E ■ B ^ 0, then it is possible to make a Lorentz transfor- 
mation to a frame in which B and E are parallel (or antiparallel) . In 
this case we can take a = E, and b = B [recall that we are using 
Heisenberg's original notation for a and b, as in (1.5)]. Then we can 
rewrite the effective Lagrangian (1.2) as 



-spinor 4^3 A^ J f n \l + ieBtY \B 

n=l V \ / \ 

' ^ ' coth ^nvrU (1.26) 



\-{eEtf \E 

using a remarkable trigonometric identity due to Ramanujan and 
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Sitaramachandrarao [27, 28] 



^2 °2^ coth f n7r| 

[ttx cot TTx) {Try coth vry) = 1 H (y^ — x^) — 27rxy^ > — — :5-^ ^ 

3 '^^ n(n^ + 7/^1 



oo 



coth (niT 



y 



-27rx^yT , ; %{ ■ (1.27) 



?i=i 



Thus, we see that in this case of parallel electric and magnetic fields, 
the imaginary part of the effective Lagrangian is 



lmL^.„-„„^ ~ „ ^ > — coth I -^ nvr I exp 



-spinor 87^2 Z^„ \^ ^ 

n=l 



m?"im 
eE 



fl.28) 



a result stated by Bunkin and Tugov [40], and obtained in a very 
different manner by Nikishov [29] using a virial representation of the 
imaginary part of C. Note that as i? ^ we recover the Schwinger 
instanton result (1.25), and that for i? > we get a slight, but not 
strong, enhancement to the prefactor, with the exponent being un- 
changed. 

1.3.2. Scalar QED case 

The analysis of the previous section for the weak field expansions of the 
spinor QED effective Lagrangian (1.2) can also be applied to Weisskopf's 
effective Lagrangian (1.8) for scalar QED. The weak field expansion of (1.8), 
expressed in terms of the Lorentz invariants a and h defined in (1.5), is: 

(1) m^ \^(2 oy Y^ B2kB2n-2k ('^^^\ f 2iea \ 

^scalar ~ ig^2 Z^i ^ ^" ^^ (2fc)!(2n - 2A:)! V "^V \m? ) ' 



n=2 A;=0 



(1.29) 



Where we have defined the convenient shorthand 



B2k = ' ^2k-i ^ 2fc , (1.30) 

in terms of which the relevant expansions are [compare with (1.16)]: 

sinhz ^ (2A;)! ' sinz ^ (2A;)! ' ^ ' ^ 

fc=o ^ ^ fc=o ^ ^ 

As in the spinor case, at any given order of (1.29), the combinations of a and 
b can be rewritten in terms of the more familiar Lorentz invariants {E'^ — B^) 



14 Gerald V. Dunne 



and E ■ B. For example, the first two orders yield: 

(1) '^ \360 36 "•" 360^ '^ \^ 2520 """ 360 360 """ 2520^ 

scalar IQm^TT^ ^ IGm^TT^ ^ ' ' ' 

_ e^[7{a^-b^)^ + 4{abf] e^a^ - b^) [31{a^ - b^^ + 77{abf] 
~ 5760m%2 "^ 40320m%2 ^ '" ' 

(1.32) 

Tlie first term in (1.32) agrees with the first term listed in [2]. 

Several special cases of the general expansion (1.29) are of interest. 

(1) Purely magnetic background : If the background is purely magnetic, 
of strength B > 0, then the integral representation (1.8) reduces to 

>c£L = ^ /"^ ^ f-V - - + ^1 ^""^^^'^^"^ (1-33) 

scalar -|^g^2 J^ g2 ysmh S S QJ 

which has the asymptotic expansion 

,(1) _ m"^ y^ ^2n+4 ^ 26^ ^""^*^ 

n=0 



scalar iQt^2 Z^ ^2n + 4) {2n + 3) {2n + 2) \ m^ J ' ^- ' 



Once again, this is a divergent but Borel summable perturbative series, 
and the Borel sum gives precisely the integral representation (1.33). 

(2) Purely electric background : If the background is purely electric, of 
strength E > 0, then the integral representation (1.8) reduces formally 
to 

(1) ^e^ j^ds_ / J_ _l_s\ 2,/(,^) 

^scalar ig^2 _/^ ^2 l^g^^S S q) ' ^ ' 

This integral representation has the asymptotic expansion 

(1) n^ ^ (-l)"^2n+4 (2eE\ 2"+^ 

scalar ~ ;Lq^2 A^(2n + 4)(2n + 3)(2n + 2) V"iV ' ^' ' 

As in the spinor case, this is a non Borel summable perturbative series, 
and correspondingly the effective Lagrangian has a nonperturbative 
imaginary part given by 



2tp2 ^ t ^\n-l 



^"i^scalar ~ ^6^ 2^ „2 «^^P 

n=l 



.2„ 



(1) e^£'^ Y^ (-1)'" " m"7rn 



eS 



(1.37) 
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(3) Perpendicular magnetic and electric background fields : li E ■ B = 0, 
then Q = 0, and from (1.6) we see that either a = or 6 = 0, depending 
on the sign of J^. If E'^ > B^, then the situation is as for the purely 
electric field, with a vacuum instability. On the other hand, if B^ > 
E"^, then the system behaves like the case of a purely magnetic field. If 
E"^ = B^, then CgJ^i^^j. vanishes and there are no quantum corrections 
at all. 

(4) Parallel magnetic and electric background fields : If E and B are par- 
allel, we can rewrite the effective Lagrangian (1.8) as 

,(1) _ e^EB ^ p dt e-™'"^* f [eBtf (-1)" 



scalar 



8^^ ^Jo t n 1 1 + (e5t)2 sinh (I 



riTT 



{eEtf (-1)" 1 

+ l-(eii;t)2sinh(|n^)J • ^^"^^^ 

using another remarkable trigonometric identity due to Ramanujan 
and Sitaramachandrarao [27] 

vrx \ / Try \ ^ , vr^ 2 2^ 



l + —{x'-y' 



sin-TTX/ \sinh7ry/ 6 

1)" 1 o 3V^ (-1)' 






^ n(n2 - x2) sinh (n7r|) ^ n(n2 + y2) ^^^^ ^ 

^1.39) 



nvr- 



Thus, we see that in this case of parallel electric and magnetic fields, 
the imaginary part of the effective Lagrangian is 



(1) e'EB ^ (-1)"^! 1 



eE 



:i.40) 



a result first obtained in a different manner by Popov [30] . Once again, 
note that as i? ^ we recover the pure electric field result (1.37), and 
that for S > we get a slight, but not strong, enhancement to the 
prefactor, with the exponent being unchanged. 
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1.4. Strong-field expansions of Heisenberg-Euler 

In order to discuss the strong-field expansions of the Heisenberg-Euler ef- 
fective Lagrangians (1.2) and (1.8), it is useful to convert the proper-time 
integral representations into zeta function expressions and alternative in- 
tegral representations which are more suited to the strong field limit. We 
concentrate here on the case of a magnetic background field. In zeta func- 
tion regularization [5,6,31,32] we define the determinant of an operator O 
in terms of its (^-function 

c(«) = E^~' (1.41) 

A 

where the sum is over the eigenvalues A of C Then the determinant of O 
can be defined as 

detC' = exp[-C'(0)] . (1.42) 

For the Heisenberg-Euler effective Lagrangians in a magnetic background, 
these ^-functions can be expressed in terms of the standard Hurwitz (- 
function of number theory [22,23,25]: 

oo 

CH{s,z) = Y,in + zr', Re{s)>l, ^ / 0,-1,-2,... . (1.43) 

ra=0 

This Hurwitz ^-function can be analytically continued in the s plane to define 
an analytic function with a single simple pole at s = 1. 

1.4.1. Spinor QED case 

For spinor QED in a magnetic background, the eigenvalues of the Dirac 
operator are 

\^=m^ + k'i + eB{2n + l±l) , n = 0, 1,... (1.44) 

where the it refers to the two spin components, and k± is the transverse 
2-momentum. Then the ^-function for this Dirac operator is [5] 

Cspinor(.) -^2^2^J J^ [ -2 ; 



n=0 ± 

m^ UbV (A) 

47r2 \w? ) {s - 1) 



2Qh\s-1, 



2 \ / 2 \ ^ — s 

m \ I m ^ 



2eB / V 2eB 



[I 



:i.45) 
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Here the overall factor of ^ is the Landau degeneracy factor, and the renor- 
malization scale n has been introduced to make the eigenvalues dimension- 
less. Given this form in terms of the Hurwitz zeta function, it is straightfor- 
ward to obtain 



~>spinorv'-^7 




(1.46) 



With on-shell renormalization (p = m), and subtracting the zero field con- 
tribution. 



■ 4^ , which ensures that C: 



spinor 



vanishes when B = 0, we obtain 



C 



(1) 



spinor 



{eBf 

27r2 



{cK-.^)-c«(-i,^).n 



1 1 



12 



+ 



2^) 



m 
2^ 



(1.47) 



Here we have used the fact [22,25] that Ch{—^,v) = — ^^^ + f ~ T' 

We can exhibit directly the equivalence of this ("-function representation 
(1.47) with the proper time integral representation (1.18) by noting the fol- 
lowing integral representations of the Hurwitz (-function [22]: 



Ch{s,z) 



1 



r(5) JO 

l-s 






1 



dt 



Re{s) > 1 , Reiz) > 



z 



z ^ sz ^ '^ 

1+— +^^ 



+ 



«-l 



T s) 



dt 
7T^' 



'^^Ucotht---- 






3, 

1.48) 



where the second expression is valid for Re{s) > —2, as a result of the 
subtractions of the integrand [34]. Thus, we can evaluate the derivative at 
s = — 1, as required: 



1 



Ch{-'^,z) = tt;- -r -CH{-'i^,z)lnz- - -^ e 



12 



dt 



-2zt 



1 t 

coth t 

t 3 



(1.49) 



Comparing with (1.47), we recover exactly the familiar proper-time integral 
representation (1.18) for the one loop Heisenberg-Euler effective Lagrangian 

(1) 

"spinor' 



c 
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Now to make a strong-field expansion of /^gp-jjoj, we use the following 
relation between the Hurwitz (^-function and the log of the T function [22,25]: 

C^(-l, z) = C'(-l) - I ln(27r) - |(1 - z) + y" lnr(x)dx . (1.50) 

This identity follows from an integration of Binet's integral representation 
[25] of Inr(z). Thus we can write 



C 



(1) _ jeBY 

spinor 27]-2 



+ 



1 , m? 3 / m^ \ m? 



ln(27r) 



1 



+ 



2 1 / 2 \ 2' 

m 1 I m ^ 



12 4eB 2 \2eB 



rn?' \ r 2eB 



(1.51) 



In the strong-field limit, the range of integration in the last term vanishes, 
so we can use the Taylor expansion [22,23] of lnr(a:): 



In r(a;) = —Inx — jx + N^ 



-ir 



n=2 



n 



■C{n)x'' 



(1.52) 



where C{n) is the usual Riemann (^-function. This leads to : 



C 



(1) ^{eBf 

spinor 27]-2 



ra" 3 / ra 
+ 



12 "^ ^ *- ^' 4e5 ' 4 V 2eB 



m 
4^ 



ln(27r) 



+ 



2 1 / m^^ '■ 
12 ' AeB 2 V 2^6 



1 m 



2Ib) ~ 2 \2^) 



m 



2eB 



1-ln 



m 



2eB 



n=2 



nin + l) \2^) 



2 \ "+1' 



(1.53) 



The leading behavior in the strong field limit is 



C 



(1) 



''^'^nr^~,+ 



-spinor 247r2 V "^^ 

which agrees with Weisskopf's original observation (1.12). 



(1.54) 
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1.4.2. Scalar QED case 

For scalar QED the analysis is very similar. There is no spin projection term 
in the eigenvalues of the Klein-Gordon operator, so the ("-function is 

, ._eB -^ r d'^k^ (m^ + k\ + eB{2n + l)" 

Cscalar(^) -^2^J J^ [ -2 

= ^(2^) (73T)C4-i>2 + 2^J- ^'-''^ 

With on-shell renormalization (p = m), and subtracting the zero field con- 
tribution ^^, we obtain 

-^1- 47r2 P^ V '2^ 2eBj^ i\2eBj 



+ 



1 + ln 



/ m 






2eB 



-.2 



c«i-i.^ + ^ir (1.56) 



We can exhibit directly the equivalence of this ("-function representation to 
the proper time integral representation (1.33) by noting the following integral 
representations of the Hurwitz ("-function: 

1 1 f°° e~*/^ 1 
Cj,{s,- + z) = -— e'^H'-'- -dt , Re{s)>l,Re{z)>-- 

2 r(s) Jo 1 - e""* 2 

z^-' sz-^-' 2*-i /""^ dt _2,t ( ^ 1 * 



1 24 r(s) io t^'" \sm\\t t 6^ 

(L57) 

where the second expression is valid for Re{s) > —2 as a result of the sub- 
tractions of the integrand. Thus, we can evaluate the derivative at s = — 1, 
as required: 

1 1 z^ z'^ 

Comparing with (1.56), we recover exactly the familiar proper-time integral 
representation (1.33) for the one loop Heisenberg-Euler effective Lagrangian 

scalar ' 
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Now to make a strong-field expansion of C^J^^i^j. we use the following re- 
lation between the Hurwitz ^-function and the log of the T function [22,25]: 

C^(-l,i + ^) = -k'(-l)-Jln(2vr)-i^ + ^+ f\nr(x + l)dx 

:i.59) 



2 ' 2 ' ' ' 2 ' ' 24 2 7o V 2 



where we have used the fact that [33,34] 



\nT{x)dx = A ln2 + ^ ln7r + ^ " ^C'(-l) • (1-60) 

The strong-field expansion is obtained using the Taylor expansion: 

lnr(z + -j =- ln7r-(7 + 21n2)a; + ^^ '- ^ ^- ((n) x"" . 

n=2 

(1.61) 

Thus, we obtain the strong-field expansion 

"^ _ 1 / m2 ^ =^1 
24 ~ 2 V2aB 



.(1) _ ieBrf5fm^\' 

scalar 4^2 ) 4 I 2eB ' 



1 + ln 



2 
m 



2eB 



,2 \2 



1 >// N In 2 m , , , 1 f m^ \ , 

-C'(-l) n(2 - - (7 + 2 In 2) 

2^ ^ ^ 24 4eS ^ ^ 2 \2eB J ^' ^ 



+E 



n(n + l) \2eB 

n=2 \ > / \ 



(1.62) 



The leading behavior in the strong field limit is 

4;L~yri"(7;^)+-- ("3) 

2. Inhomogeneous backgrounds : beyond the constant field 
strength limit 

The Heisenberg-Euler and Weisskopf results, (1.2) and (1.8) respectively, 
are computed for the special case where the background field strength F^^^ 
is constant. Clearly, it is of interest to be able to compute the effective 
Lagrangian in more general cases, where F^^ is inhomogeneous. The general 
formalism for computing such determinants was developed by Schwinger 
[12, 13], Salam and Matthews [35], and Salam and Strathdee [36], among 
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others. It is, unfortunately, cumbersome to do such calculations, and so 
several physically relevant approximate techniques have been developed. I 
briefly review these in this section. 

2.1. Solvable inhomogeneous backgrounds 

There exist a limited number of solvable cases in which one can compute 
the effective Lagrangian in a more-or-less closed-form, analogous to the 
Heisenberg-Euler result (1.2). These solvable cases reveal some interesting 
new physics, beyond the constant field approximation, as well as provid- 
ing useful checks on the more general approximate techniques. In order 
for a particular background to lead to such an explicit form for the effec- 
tive Lagrangian, we need to be able to solve the spectral problem for the 
corresponding Dirac or Klein-Gordon operator appearing in (1.1) or (1.7), 
respectively. Solvable Dirac and Klein-Gordon equations have been stud- 
ied exhaustively, so one can simply look through the list of solvable cases 
and decide if the corresponding electromagnetic background is physically 
relevant. 

2.1.1. Plane-wave background 

The simplest case of an inhomogeneous background is that of a single 
monochromatic plane wave field. This was shown to lead to a solvable Dirac 
equation by Volkov in 1935 [37], and is described in [38]. Such a gauge field 
can be written A^ = A^{k-x), where k^ is light-like (k-k = 0), and k-A = 0. 
Then the solutions to the Dirac equation can be written explicitly as 



V'p = e' 



in 



e(7-fc)(7-^) 
2{k-p) 



^(^) (2.1) 



^2po 



where p is the momentum, u{p) is a fundamental free spinor, and the phase 
factor O is 

Having the exact solutions means one can compute the exact Green's func- 
tion, which in turn means one can compute the one-loop effective action. 
Schwinger [12] showed that the effective action in fact vanishes for such a 
plane wave background. (This is consistent with the fact that for a plane 
wave both Lorentz invariants J-' and Q vanish, so that the low frequency 
limit Heisenberg-Euler result vanishes.) In Schwinger's words [12]: "there 
are no nonlinear vacuum phenomena for a single plane wave, of arbitrary 
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strength and spectral composition". In particular, Schwinger's result means 
that there is no pair production in a single monochromatic plane wave back- 
ground, which can be understood from simple energy and momentum con- 
servation - at least two plane waves are necessary to produce a pair [12,39]. 
The problem of pair production in crossed laser beams has been addressed 
recently in [41]. 

2.1.2. Time dependent electric field E{t) = E'sech {ujt) 

With such an electric field, pointing in a fixed direction, spatially uniform, 
but with a special prescribed sech {ujt) time dependence, the Dirac and 
Klein-Gordon equations reduce to hyper geometric equations. Narozhnyi and 
Nikishov [42] (see also [30]) used this fact to express the effective action as 



ImS = ± \, .„^ / dV / dp, log (1 ± wp) (2.3) 



(2^ + 1) I ,2. 
2(27r)3 

where the upper /lower signs it refers to spinor/scalar QED, respectively, and 
(2s -|- 1) is an overall spin factor of 2 for spinor QED and 1 for scalar QED. 
The momentum has been divided into pz along the direction of the electric 
field, with pj_ being the transverse momentum. The Wp can be computed as 
a scattering amplitude factor for the corresponding Dirac or Klein-Gordon 
equation. For spinor QED one finds [42]: 

spinor siuh 7r(A — ^ + u) sinh 7r(A + ^ — v) 



^ sinh7r(— A + ^ + z/) sinh7r(A + /[/ + jy) 

where A = — ^, and 



(2.4) 



^^ = —\lm?+p\_ + {pz + eE/u:) 



— ^m2 +p^ + {pz - eE/ujf . (2.5) 



For scalar QED one finds [42] 



^scalar ^ cosh 7r(A^ - II + i^) cosh 7r(A^ + H-v) ,^ g. 

^ cosh 7r(A' — fj, — u) cosh 7r(A' + /j, + i') 



where A' = y^A^ — 1/4, and fj, and ly are as defined in (2.5). 

When (J ^ 0, the time scale of the electric field tends to infinity and we 
recover the constant field results (1.25) and (1.37). To see this, note that 



Vp 



l=F??p 



vr 



where ??„ = exp ^("t- +P±) ■ (2-7) 

eE 
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(±1)" 



rjp from the expansion 



We get a factor of eE from J dpz, a factor of 

of the log in (2.3), and finahy a Landau degeneracy factor of ^ from J (Pp±. 

Thus, the virial representation (2.3) yields: 



2j?2 



ImC = (2s + 1 



167r3 



E 

n=l 



(±1) 



n— 1 



n^ 



exp 



m?TTn 
eE 



(2.8) 



in agreement with (1.25) and (1.37). 

For nonzero frequency uj, the Narozhni-Nikishov answer (2.3) is expressed 
as the integral representation (2.3) with Wp given by (2.4) or (2.6) for spinor 
or scalar QED respectively. In [43], for spinor QED, a more explicit form of 
this exact answer was found, involving just a single integral representation, 
analogous to (1.23) for the constant field case. 



S 



"3^^2— 



dt 



UJ 



eE tuj' 



o2nt 



where f ^ 



fu)^ 



m^ 



')2F^(l 



1 3 



2' 2' 1 



+ {E 



-E) 



(2.9) 



+ 2^t. This can be compared with the proper-time repre- 



sentation of the constant field result (1.23), which can be rewritten as 
(e^)2 f'^ ds 



C 







AeE 



■t 



m^ 



+2 arctanh 



(2.10) 



2.1.3. Space dependent magnetic field B{x) = i?sech (^) 

As in the previous case, this case is soluble because the Dirac and Klein- 
Gordon equations reduce to hypergeometric equations. In [44] an exact 
closed-form expression was found for the spinor effective Lagrangian in 2 + 1 
dimensional QED in such a background, and in [45] this was extended to 
3+1 dimensions. In the 2 + 1 dimensional case [44], 



s: 



(2+1) 



L 

47rA2 



dt 



o2nt 



±1 



(6± - it) 



(A^m^ + vi), Xm-tv± 
In- : — 

v± Am + iv± 



+ c.c. 



(2.11) 
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where + denotes scalar QED, — denotes spinor QED, c.c. denotes the com- 
plex conjugate, and 

V(eBA2)2 + 1/4 (+) bosons 
eBX^ (— ) fermions 

^2 ^{t^ + 2ith+- 1/4 (+) bosons 

^^ \t^ + 2ith^ (-) fermions ^ ' ' 

For spinor QED this was extended to 3 + 1 dimensions [45] , to give 



b^ = { v^;:r^ ^ ' vi ::rr^^ (2-12) 



^(3+1) 2m^L^\ r dt 



spinor 



3vr2 Jo e^''* - 1 



3 



(eB it \ {l + v^)i 
— o oTo I arcsm(zt;) + c.c. 



(2.14) 



where v"^ = — ttt + 2i^ t. This can be compared directly with the constant 
magnetic field Heisenberg-Euler result (1.18), which can be written in a 
modified form: 



(eBf f^ ds _^J . 1 

- „ ' / -^e eB'' coth s 

8vr^ 7o s V s 


~3J 

\ 

i In 


1 + i 


f2eBt\ 


+ C.C. \ . 

(2.15) 


87r2 \m^ ) Jo e'^^t - 1 i ^v "i^ 



2.1.4. Other solvable backgrounds 

Another solvable case, which is very similar to the cases described in Sections 
2.1.2 and 2.1.3, is that of an electric field pointing in a fixed direction, static 
in time, but having spatial dependence 

-E;(x) =^sech2(^) (2.16) 

A 

along the direction of the electric field. Once again, this case is solvable 
because with the gauge field ^0(2;) = —EX tanh(^), the Dirac and Klein- 
Gordon equations reduce to hyper geometric equations. For details see [46] 
and [47]. Such a background in Euclidean 2-dimensional space was studied 
in [48]. 

An interesting class of solvable electric fields was recently found by 
Tomaras et al [49,50], who found that when the electric field depends on 
a light-cone coordinate x^ it x^, then the one loop effective Lagrangian can 
be found exactly, and has essentially the same form as the constant field 
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Heisenberg-Euler result. Specifically, define the light-cone coordinate 

x+ = ^{x^-x^) (2.17) 

and choose the vector potential 

A^{x+) = / duE{u) (2.18) 

Jo 

with A± = 0, and use the gauge A+ = 0. Then the spinor one loop effective 
action is [49,50] 

CL = ^2J<^"^ [% ^"^^™' [^E{x^) « -th {eE{x^)s) - l] . 

(2.19) 

This leads to a pair production rate 

e'^E'^{x+) 1 ^ feEix^ 

n=l 

(2.20) 



'P y I 247r 47r ^ 



exp 



mrm? 



|e^(x^ 



2.2. Semiclassical approximation 

The examples listed in the previous section exhaust the analytically solvable 
cases with physically realistic inhomogeneous electromagnetic backgrounds. 
For more general forms of inhomogeneity in F^j^i, a powerful approach is the 
semiclassical approach. The semiclassical method works best if the back- 
ground field depends on just one coordinate. Here I concentrate on the 
example provided by an electric background field which points in a fixed 
direction in space (say the z direction), but has a time dependent profile. 
This approach was pioneered by Brezin and Itzykson [51], and Popov and 
Marinov [30,52], building on earlier ionization work by Keldysh [53] and 
Perelomov et al [54] which will be discussed below. 

For a time dependent electric field we can choose the gauge 

A,{t) = -- f{iot) ^ E,{t) = EJioot) . (2.21) 

For example, the solvable case in Section 2.1.2 corresponds to the choice: 
■Az{t) = —— tanh(a;t). I also concentrate in this section on scalar QED; 
spinor QED adds some technical details concerning the spin components. 
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but the basic techniques and results are very similar. Then the Klein-Gordon 
operator for such a background field is 



We expand the field operator in solutions of [51] 

-4> - {p, - eA,f(t) = (m2 +p\)(l) 
with scattering boundary conditions 



(2.22) 



(2.23) 



~ e"**"*^™^^^^ + b^e^^^™'^^'^^ 


t - 
t- 


-> — oo 


~ ap^e-^*V-^+P' , 


-^ +00 



(2.24) 



so that particles are viewed as propagating forwards in time and antiparticles 
as propagating backwards in time. Then the pair-creation probability is 
given in terms of the refiection coefficient b^ 



dPp 



3 I^PI 



(2.25) 



This is an example of "over-the-barrier" scattering [55] , and the semiclassical 
WKB expression for this amplitude gives 



exp 



-2 Jm * Jm"^ +p]_ + \Pz- eA2(t)]2 dt 



(2.26) 



Changing variables to u = i[pz — eAz{t)]/m, we can write 



exp 



m 



1 — V? + p^^/m? 



\E{u) 



du 



(2.27) 



where the limits of the u integral are set by the turning points. The ex- 
ponential term for the pair production rate comes from neglecting the mo- 
menta, which instead affect the prefactor (in fact, different prescriptions 
for treating the momenta lead to slightly different forms of the WKB pref- 
actor [30,51,52].) The leading WKB approximation for the pair-creation 
probability (which is directly related to the imaginary part of the effective 
Lagrangian) is then 



P~ exp[-7r^c5(7)] 



(2.28) 
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where Ec = 
Euler, and 



eh 



is the critical electric field (1.3) defined by Heisenberg and 



2 /"I ^i - n2 
vr /_i LB u 



Here 7 is the " adiabaticity" parameter : 



7 = 

For example, in the sinusoidal case, 

E 



eE 



^. 



UJ 



sin(c<;t) 



E^{t) = Ecos{ujt) 



(2.29) 



(2.30) 



(2.31) 



we can evaluate (7(7) exactly [30,51,52], 



Egi-f) = - 




E 



7 



1 + 7^ 



(2.32) 



where K(x) and E(x) are the complete elliptic integral functions [25]. 

This is a truly remarkable result as, when inserted into the pair produc- 
tion rate (2.28), it leads to an interpolation of the pair production probabil- 
ity between a nonperturbative limit when the adiabaticity parameter 7^1, 
and a perturbative limit when 7^1: 



P~ < 



exp 



-TT 



eEh 



g^ \4mc^/hijj 



V. \u}mc/ 



7^1 (nonperturbative) 



7^1 (perturbative) 



(2.33) 



where we have re-inserted the h and c factors. When 7^1 the time 
dependent field is slowly varying and we recover the leading nonperturbative 
exponential behavior of the Heisenberg-Euler result (1.10). On the other 
hand, when 7^1 the result in (2.33) has the perturbative form of the 
square of the gauge field strength (^)^ raised to a power equal to the 
number of factors of the photon energies fko needed to make up the pair 
creation threshold energy 2m(?'. 
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It is worth pausing here to recah the corresponding nonrelativistic ion- 
ization problem, where this type of interpolation behavior was first found 
in a classic paper by Keldysh [53]. Oppenheimer [56] first estimated the 
ionization probability for a hydrogen atom in a static electric field E to have 
the exponential form: 



P 



exp 



3 Eh^ 



(2.34) 



More generally, for tunneling through a barrier with binding energy £i„ the 
WKB expression is 



P 



exp 



exp 



2 

'h 

4 
'3' 



St/eE 



dz\/1m{E\i — eEz) 



^2mS, 



3/2 



eEh 



(2.35) 



For a hydrogen atom, S^ = ^^, leading to Oppenheimer's estimate (2.34). 
In a seminal paper [53], Keldysh generalized this ionization analysis to a 
time dependent electric field, and in particular a monochromatic sinusoidal 
field E{t) = Ecos{ujt). This approach was further developed by Perelomov 



/2m£, 



3/2 



et al in [54]. With the weak field condition, E <^ — ^ — , and the condition 
for a classical electromagnetic field, huj <^ S^, Keldysh defined the ratio of 
these two small ratios to be the " adiabaticity parameter" : 



7 



hui/Sb 



ehE/{y/2^£ 



3/2 ^ 



ijj\J2m£i) 
eE 



(2.36) 



Using WKB he computed the leading behavior of the ionization probability 
to be 



P 



exp 



2^ 



5(7) 



(2.37) 



where for the sinusoidal electric field 

1 



ail) 



1 + 



77 I arcsinh7 
27 V 



27 



[§7, 
\ln(27) 



7« 1 
7> 1 



(2.38) 



Keldysh noted that the ionization probability (2.37) interpolates, as 7 ranges 
from 7 <C 1 to 7 ^ 1, between the nonperturbative tunneling form (2.35) 
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and a perturbative multiphoton ionization form: 



P~ < 



exp 



4 V2m£f 



3/2 



f eE 



eEh 



2ef,/hjj 



^ \ 2uiyJ2rr£\y 



7 <C 1 (nonperturbative) . 



7^1 (perturbative) 



(2.39) 



For 7 <C 1 we recover the tunneling form (2.35), and for 7 ^ 1 we see a 
multiphoton form of the ionization rate. The analogy of Keldysh's result 
(2.39) to the QED case for the WKB vacuum pair-production result (2.33) 
becomes clear taking the "binding energy" for the pair production process 
to be the pair rest energy : Ei, = 2ra(? . 



2.3. Derivative expansion 

Another approach to the problem of an inhomogeneous background field 
strength F^^ is to expand about the solvable constant field case, in terms 
of derivatives of F^i,. Thus, generically one obtains an expression for the 
one-loop effective action of the form 

S = 5[o] [F] + S[2] [F, {dFf] + 5[4] [F, {dFf, {dF)^] + ... (2.40) 

where S\2n\ involves up to 2n derivatives of the field strength. This formal 
expansion should of course be understood as a multi-dimensional expansion, 
because even the constant field term S'ro][-F], which is the Heisenberg-Euler 
result, has itself a doubly-infinite weak-field expansion, as shown in (1.15). 
At higher orders in the derivative expansion, things rapidly become more 
complicated since there are many more Lorentz and gauge invariant terms 
that can (and do) appear once we include derivatives of F. In the next 
two subsections 2.3.1 and 2.3.2 we consider the question of the convergence 
of the derivative expansion and show that it should be understood as an 
asymptotic expansion. 

Nevertheless, the derivative expansion approach is a very powerful one. 
The leading order approximation consists of taking the Heisenberg-Euler 
constant field result for the effective Lagrangian, substituting the inhomo- 
geneous F^y for the constant one, and then performing the spacetime inte- 
grations. In many cases, this gives an excellent approximation to the full 
effective action. We expect the derivative expansion to be "good" (in the 
sense of an asymptotic expansion [57]) when the scale of variation of the 
background field is large compared to the electron Compton wavelength, or 
when the scale of variation of the background field is large compared to the 
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length scale set by the average magnitude of the background field. For ex- 
ample, for a spatially varying magnetic background, with variation scale A, 

and peak field B (which defines a "magnetic length" scale Ib = \/ ^)) the 
derivative expansion is expected to be good when 



A > /b or A > Acompton = — ■ (2.41) 

mc 

A very useful tool for systematizing the derivative expansion is to use the 
Fock-Schwinger gauge choice 

x^,A^ = . (2.42) 

In this gauge it is possible [19, 20] to express the gauge field in terms of the 
field strength and its derivatives (evaluated at some reference point x = 0): 



1 

(n. - 
n=0 



[n + l)n\ 



= x" / d7?7?e''^-^F„^. (2.43) 

JO 

In this gauge one can expand the effective action (1.1) in powers of deriva- 
tives of the field strength. A very convenient formalism for doing such an 
expansion is the world-line formalism [58-62] in which the effective action is 
written as a quantum mechanical path integral. For scalar QED, 



S,,[A] = I'^^e-"''^ JvxexpLj dr 



-x^ + ieA^{x{T))±^{T) 



(2.44) 



where the path integral is over closed spacetime loops x^(r). Keeping just 
the first term in the derivative expansion (2.43) of A^, this quantum mechan- 
ical path integral is Gaussian and exactly solvable. Higher order terms in 
the derivative expansion are then generated by expanding the non-Gaussian 
interaction terms and computing the resulting expansion terms using the 
known worldline Green's functions for the constant background field about 
which one is perturbing. For spinor QED one must also include Grassmann 
valued fields V'^(t) which are antiperiodic in the propertime parameter r: 
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V'm(O) = — V'm(^)- Then for spinor QED 



5sp[A]=-^r^e-^^/Px |V^l; 



^ JO 
exp { — dr 



1-2 1 

-x^ + - 

4 2 



7X^ + -^p■1p + ieA^x^" - ieipf" F^.^ip'' 



(2.45) 



which can once again be expanded systematically about the constant field 
case. The resulting next-to-leading order (NLO) derivative expansion terms 
for 2 + 1 dimensions can be found in [63], and for 3 + 1 dimensions in [64,65]. 
The general expressions, even at the next-to-leading order level, are too 
long to be usefully included here. But for a given inhomogeneous back- 
ground it is straightforward (albeit tedious) to evaluate the first few orders 
of the derivative expansion. Perhaps a more important question is to ask : 
is the derivative expansion convergent? As it stands, this is not a well-posed 
question because of the proliferation of terms with completely new Lorentz 
structures at higher orders. However, this question can be made well-posed 
if we restrict the background field to depend on only one spacetime coordi- 
nate, with the inhomogeneity characterized by a single scale. In this case 
the general derivative expansion expression (2.40) reduces to a double sum, 
which we can hope to analyze in more detail. To be even more specific, 
consider some of the solvable inhomogeneous cases from Section 2.1 - by 
comparing the closed-form answers with the derivative expansion we learn 
something about how the derivative expansion behaves. 



2.3.1. Spatially inhomogeneous magnetic field B{x) = Bsech (^) 

For this spatially inhomogeneous magnetic background field the exact closed 
form answer for the spinor QED effective action was given in (2.14). It is a 
simple exercise to make a double asymptotic expansion of this answer: 

8^3/2 Z> J! (^A)2i ^ r(2A; + l)T{2k + j + \) W^ 

(2.46) 

(where the j = and k = 1 term is excluded from the sum as it corresponds 
to a charge renormalization term.) The sum over j corresponds to the deriva- 
tive expansion, with expansion parameter -, — jrj = ( compton -j2^ while the sum 
over k corresponds to a perturbative expansion with expansion parameter 
■^, as in the constant-field Heisenberg-Euler case (1.19). The j = term 
in (2.46) is the leading order derivative expansion term, and the j = 1 term 
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is the NLO term. The LO effective Lagrangian is just the Heisenberg-Euler 
expression (1.18) for a constant magnetic field, 



m^ ^ B2k f2eB^^'' 



87r2 ^ 2k{2k - l)(2/c - 2) \ m"^ 



(2.47) 



while the NLO effective Lagrangian for a static but spatially inhomogeneous 
magnetic field is [63,65] 



s 



2 djBdjB y. B2k+2 (2eB\ ^'"'^ 
4^2^2 2^ (2fc - 1) \m^ J 



e'-n^y T^l''^ -^ (2-48) 



which involves two derivatives of the magnetic field. To compare with the 
expansion (2.46) of the exact result, we substitute the inhomogeneous mag- 
netic field, B{x) = i?sech (^), for B in the derivative expansion expressions 
(2.47) and (2.48), and then integrate over the inhomogeneity direction x. 
This yields the LO and NLO derivative expansion results for the effective 
action as 

L'^Xm^^ 1 B2kr{2k-2) f2eBY'' , ^ 

^^'^--^^l^Yk r(2fc + i) \-^) ' ^'-''^ 

LW ^ g2fc+2r(2fc-l) (2eB\^'' 

^[^i^-^A^^ r(2fe + i) \-^) ■ ^'-'"^ 

These results agree exactly with the j = and j = 1 terms from the asymp- 
totic expansion (2.46) of the exact result. This example clearly illustrates the 
approach of substituting the inhomogeneous fields for the fields appearing 
in the derivative expansion expressions for the effective Lagrangian. 

2.3.2. Time dependent electric field E{t) = i^sech (cjt) 

Another interesting comparison of the derivative expansion can be made for 
this time dependent solvable case, whose exact solution was given in Section 
2.1.2. One can do the same thing described above for the inhomogeneous 
magnetic case, and compare the asymptotic expansion of the exact result 
(2.9) with the LO and NLO field theoretic derivative expansion results from 
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[65]. Instead, here I show [66] how the derivative expansion agrees with the 
semiclassical WKB results of Section 2.2. 

An asymptotic expansion of the exact result (2.9) yields the derivative 
expansion 



S 



m 



j k ^ ' ' 



f-f- (-ly+fc /2eij;y" r(2fc + j)r(2fc + i-2)g2fc+2i 

87r3/2u; ^ 4^ (m/u;)2i \ m? ) 



j\{2k)\T{2k + j + 



(2.51) 



Fixing the order j of the derivative expansion, the remaining sum over k is 
divergent and nonalternating, and hence has a nonperturbative imaginary 
part. For fixed j, the expansion coefficients behave for large k as 



Si) 



-i)^+fer(2fc + j)r(2fc + J + 2)g2fc+2j+2 

r(2fc + 3)r(2A; + j + |) 



r(2A; + 3j 



A; ^ oo . 



(2vr)2i+2fe+2 

Using standard Borel dispersion relations [67,68], one finds [66] 

1 /vrm^u;^^"' 



ImS'^^^ 



exp 






Remarkably, the sum over j exponentiates to give the leading term: 

1 /muj\ 2^ 



ZmS ~ exp 



m?iT 



1 



(2.52) 



(2.53) 



(2.54) 



On the other hand, if we first considered a fixed order of the k summation, 
then the remaining j sum is divergent and nonalternating. For large j, the 
coefficients behave as 

Jfc) _ (_^\i+k 



T{j + 2k)T{j + 2k- 2)B2kA 



-2i 



Tij + i)r(i + 2A: + 1; 



22 



S-2fc 



r(2j + ik 



(27r)2i+2fc ' 
This leads to the nonperturbative imaginary part: 



J ^ oo 



ImS^''') 



2k 



1 /27reE\- , , , 

M!W J exp[-2vrW-] 



(2.55) 



(2.56) 



The remaining k sum also exponentiates, yielding the leading behavior 

Jm5~exp[-^^fl-^')l . (2.57) 
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How do we reconcile these two different expressions (2.54) and (2.57) for 
the imaginary part of the effective action? The answer can be found in the 
WKB analysis of Section 2.2 and the Keldysh adiabaticity parameter, which 
in this case is 

muj uj/m 

the ratio of the two expansion parameters in (2.51). 

For this particular time dependent electric field E{t) = E'sech (cot), we 
can return to the WKB analysis of Section 2.2 and evaluate [30] the exponent 
(7(7) in (2.29) which appears in the leading WKB pair production rate (2.28): 



2 f' VT 



u 



2 



IT J_i 1 + Tu^ 



v2 



2 |l-^ + ... , 7<1- 



i + yr+^ Mi-i + ... , 7»i 



(2.59) 



Then the WKB approximation (2.28) for the pair-production probability, 
P ~ exp[— 7rm^(/(7)/(ei?)], gives precisely (2.54) in the 7 <C 1 limit, and 
(2.57) in the 7^1 limit. In terms of the Borel resummation, it is a matter 
of competing exponential factors, with the dominant exponential being de- 
termined by the size of the adiabaticity parameter 7. Thus, this derivative 
expansion example is completely consistent with the WKB analysis of the 
imaginary part of the effective action. It also demonstrates clearly that the 
derivative expansion is an asymptotic expansion, rather than a convergent 
expansion, because if it were convergent there would be no imaginary parts 
generated from the derivative expansion. 

2.4. Large mass expansion 

The large mass expansion is very closely related to the derivative expansion, 
but is a different way of organizing the expansion. For example, consider 
the all-orders derivative expansion double series (2.46). This could clearly 
be written as an expansion in inverse powers of the mass m, with another 
series that involves different powers of the magnetic field strength and its 
derivatives: 

L'Xm^ ^ T{ini-2) 1 ^ (2ei?)^fe B2i-2k ., fi^ 

8^3/2 Z^ r(/ + i) m^i^ A2'-4fc (2A;)!(/-2fe)! ' ^ ' ' 

1=3 ^ k=-\ 
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The general structure of the large mass expansion of the effective action is 

nn= I f tI;^ /^-E^ojF] (2.61) 







T (47rr)'i/2 



n=l 



n! 



where the functions On[F] have mass dimension 2n and are expressed in 
terms of the field strength F^j^ and its derivatives. There are many dif- 
ferent ways to derive such an inverse mass expansion. Perhaps the most 
efficient [69] is the string-inspired worldline approach in which the effective 
action is expressed as a quantum mechanical path integral over a proper- 
time interval r G [0, T], as shown for scalar QED in (2.44) and for spinor 
QED in (2.45). These worldline forms of the effective action clearly separate 
the mass dependence from the background field dependence, and with the 
Fock-Schwinger gauge (2.43) lead to manifestly gauge invariant expressions 
for the functions ©^[F]. A key step in this process is finding a minimal basis 
for the On \F] , since these can clearly be rewritten in various ways using the 
Bianchi identity, the antisymmetry of F^jy, and integrations-by-parts. Such 
a minimal basis was constructed by U. Miiller in [70], and is described in 
detail in [69]. The expressions for the high order terms C'n[^] are very com- 
plicated and are typically generated using a symbolic manipulation program 
such as FORM [71]. Here I list, from [69], the first few terms C'n[-^] for the 
case of charged scalar particles in an inhomogeneous nonabelian background 
field : 

Oi = 0, 

02= -irF^xFxn , 

D 

1 2 

O3 = — ti F^x^F^^x — — i tr F^xFx^F^K , 

^4: ^ ^TT 'J^ ^ kX^ Xi^^ fMU^UK + 7^ ^^ ^ kX^v^ Xnvfi 

2 4 

~r "^rz 'J^ J^ nX't' Xk^ ^u^ufi 1 TTP' ^^ ^ kX^ XpL^ kv^ u^ 

fi 8 

17Z ^ ^^ ^ kX^ fiXu^ fj.UK TnF ^^^ KX^Xfiv^KUfi 

+ ^ tr F^xFf^uFxnFu^ ■ (2.62) 

Here, Ff^i^p = D^F^p, so it is clear that O^i-^] has mass dimension 2n. 
The next two terms, C'5[-F] and 06[-^]) are also presented in [69], and the 
dependence on a background scalar field is included also. The number of 
terms appearing in a given ©^[F] grows rapidly with n. For example, for 
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a gauge background, using the minimal basis, the number of independent 
terms in ©^[-F] is found [69] to be equal to 0, 1, 2, 7, 36, 300, ..., for 
n = 1, 2, . . . , 6, which looks like factorial growth. 

2.5. Worldline loops 

A recent interesting proposal for evaluating effective actions is to implement 
the worldline expressions (2.44) and (2.45) numerically, by doing the quan- 
tum mechanical path integral as a Monte-Carlo integration over a randomly 
generated set of worldline loops [73,74]. In principle this is a very powerful 
and general technique, since it is not constrained by the "slowly varying" 
restriction for the applicability of the derivative expansion or the large mass 
condition for the inverse mass expansion. Moreover, it can be applied to 
both scalar and spinor, and to both abelian and nonabelian systems. So far, 
the method has been tested on some fairly simple abelian models (as well 
as to the Casimir effect for surfaces of nontrivial shape [75]), and compares 
favorably with other methods when a comparison is possible, and produces 
robust answers in other inaccessible cases. Technical aspects of this method 
are still being developed, but this is an approach whose generality makes it 
very promising. 

2.6. Bounds 

Another approach to the evaluation of one-loop effective actions in nontrivial 
inhomogeneous backgrounds is to seek rigorous mathematical bounds using 
operator theory and functional analysis. This approach has been developed 
recently by M. Fry [76], and also has great potential to yield important infor- 
mation about the behavior of fermion determinants in general backgrounds. 
For example, in Euclidean 2-dimensions, for a unidirectional abelian field 
strength F{x), the fermion determinant is bounded below by [76] 

In det( -ilp + m) > — d'^x eFix) - irn^ + eF(x)) In ( 1 -^ lAp. ) 
Air J I ^ \ m^ J 

(2.63) 

For higher dimensional systems, such general bounds are considerably more 
difficult to establish, but once again their appeal is their generality. 

3. Heisenberg-Euler beyond QED 

In 1965, Vanyashin and Terent'ev [77] performed a Heisenberg-Euler com- 
putation beyond QED, in a nonrenormalizable theory of vector electrody- 
namics. They computed the Heisenberg-Euler effective Lagrangian for a 
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massive charged vector boson of gyromagnetic ratio 2, in a constant elec- 
tromagnetic background field. For an electric background field they found 
a scalar pair production rate of the form in (2.8) for spin 1. More sig- 
nificantly, they noted the anti-screening effect of the vector bosons on the 
renormalized electric charge. In hindsight, this is a precursor to asymptotic 
freedom, but at the time the physical significance was not appreciated, in 
large part due to the fact that this was a nonrenormalizable theory. Later, 
after the discovery of asymptotic freedom [78,79], Skalozub [80] explained 
the physical significance of Vanyashin and Terent'ev's results by doing a re- 
lated Heisenberg-Euler effective Lagrangian computation in a renormalizable 
SU{2) theory with spontaneous symmetry breaking. Once the principles of 
consistently quantizing nonabelian gauge theories were established, a num- 
ber of Heisenberg-Euler applications and calculations followed. 



3.1. Covariantly constant nonabelian backgrounds 

The natural nonabelian generalization of Heisenberg and Euler's condition 
of constant abelian field strength is that of a covariantly constant nonabelian 
field strength: 

D^,F,p = . (3.1) 

This condition implies that the corresponding nonabelian gauge field A^ can 
be expressed as a quasi-abelian gauge field, up to a gauge transformation 

[81-83]: 

^M = -^F^^x'^in-T'^) -^d^UU'' . (3.2) 

^ 9 

Thus we can take the field strength to point in a particular direction n"" (with 
n'^n'^ = 1) in color space, and characterized by an "abelian" field strength 
F ■ 



C = ^M- «" > ^M = - o F^^-^" ^" ■ (3-3) 



Brown and Duff [81] presented the general formalism for computing one- 
loop effective actions in nonabelian theories, generalizing the approaches of 
Heisenberg and Euler, and Schwinger. Then Duff and Ramon-Medrano [84] 
applied this to the one-loop effective action in Yang-Mills theory, finding an 
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ultraviolet finite expression 
1 „ [^ (h 



C 



YM 



2(47r) 



rTr 







g-XSg-F(y;s) 



^F{y;s) 



12 



where 



F{Y;s) = -trln 



sin(7s^ 
7s 



-X' 



1 r) yfivyfjiv 



7/11/ — ^/la^i^a 



-Y Y 



fll/^ flU 



(3.4) 



(3.5) 



and X, Y and 3^ are expressed in terms of the covariantly constant back- 
ground field strength F^^ as [84]: 



rah 



xz = 'igr^F' 



flU 



abc jpc 



Y. 



ab 
jab 



gjacbp. 



^u ^ccP 1 



acb 



V"'° = n f°-'^° F' 



flU 



flU 



(3.6) 



Duff and Ramon-Medrano noted that the one-loop effective action (3.4) is 
UV finite but infrared divergent, due to the massless nature of Yang-Mills 
theory. 

Savvidy et al [82, 85, 86] considered the case of SU{2) Yang-Mills the- 
ory and wrote a formal unrenormalized expression for the one-loop effective 
Lagrangian, motivated by the original Heisenberg-Euler form (1.2): 

1 1"^ ds gfis gf2S 



C 



SU{2) 



87r2 



smh{gfis) sm{gf2s) 



+ 



1 

4^ 



ds 



{9hs){gf2s) 







sin(5'/is) sin(5r/2s) 



sinh(5r/2s) smh{gfis) 



(3.7) 



where /i and /2 are essentially Heisenberg and Euler's a and b parameters 
(1.6), with J^ = l-F^jy-F^jy, and Q = \F^^F^j^. In the chromomagnetic case 
where fi = H and /2 = 0, the effective Lagrangian has the strong field limit 



C 



11 

'487r2 



(gHf 



In 



gH 

^2 



(3.8) 



Identifying the coefficient of the log term with the /3-function, in the spirit 
of Weisskopf's original observation [2, 15] and the work of Coleman and 
Weinberg [87] and Ritus [88-90], Savvidy et al noted that this Heisenberg- 
Euler effective action computation gives the correct one-loop Yang-Mills f3- 
function [78,79]. The negative sign of the /3-function, which signals asymp- 
totic freedom in Yang-Mills theories, arises physically because of the mag- 
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netic moment coupling of the charged vector nonabehan field. The signif- 
icance of this magnetic moment coupling had been emphasized previously 
by Tsai et al [91], as noted by Salam and Strathdee [36] in their paper 
which also helped to introduce zeta function techniques to the computation 
of Heisenberg-Euler effective Lagrangians. Savvidy et al also noted the in- 
teresting fact that for a covariantly constant chromoelectric field there is a 
vacuum instability, very much like that found by Heisenberg and Euler for 
QED in a constant electric field. 

N. K. Nielsen and P. Olesen [92] found that the situation is even more 
interesting than this, as they showed that in the covariantly constant non- 
abehan case there is in fact also an instability in a chromomagnetic back- 
ground, not just in a chromoelectric background. This is in distinct contrast 
to the QED case where a constant magnetic background has no vacuum in- 
stability. Nielsen and Olesen noted that the difference is directly due to the 
anomalous magnetic moment coupling in the Yang-Mills equations, which 
contributes to the effective action: 



S = - [dkf2 E j29H{n+'^-Ss) + k^ . (3.9) 

The mode with n = and spin aligned along the chromomagnetic back- 
ground {S3 = +1) gives an imaginary contribution to the effective action, 
indicating an instability. They further showed that this result can be estab- 
lished by a precise treatment of the physical integration contours required 
in the formal bare expression (3.7) of Savvidy et al, thereby making di- 
rect connection with the proper-time results of Heisenberg and Euler, and 
Schwinger. Subsequently, Olesen and collaborators [93] developed a color 
magnetic condensate picture of the QCD vacuum based on this chromomag- 
netic instability. This chromomagnetic instability was also pointed out by 
Yildiz et al [94], who identified an imaginary contribution in the associated 
Heisenberg-Euler effective Lagrangian, and who also computed the massive 
spinor and scalar contributions to the effective Lagrangian in a covariantly 
constant background. Then, in a development that will be of interest later in 
this review, Leutwyler [95] showed that for a covariantly constant nonabelian 
background there is an instability unless the background field is self-dual (or 
anti self-dual): F^^ = ±F^i,. At about the same time, the significance of self- 
dual backgrounds was becoming appreciated in gauge theory [96,97], having 
already been realized to be important in gravity and supergravity [98,99]. 
Dittrich and Renter [100] and Elizalde [101] further developed the zeta func- 
tion approach, and applied it to the computation of the QCD effective La- 
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grangian in covariantly constant nonabelian backgrounds. Other important 
developments have been the derivative expansion and large mass expansion 
for nonabelian theories [19,20,69,102], and also for nonabelian supersym- 
metric theories [103]. 



3.2. Super symmetric Heisenberg-Euler effective Lagrangians 

The Heisenberg-Euler approach is of great use in studies of supersymmetry 
and supersymmetry breaking. It is possible to apply the Heisenberg-Euler 
techniques, as generalized by Schwinger, Tsai and others [104] to higher 
spins, to study effective Lagrangians of super symmetric theories in constant 
background fields. This is now a vast branch of quantum field theory [105], so 
it is only possible to give a small taste of the approach here. An important 
paper in this development was by Fradkin and Tseytlin [106], who used 
dimensional reduction from d = 10 to d = 4 to compute the one loop effective 
Lagrangian for 4 dimensional AA = 4 SYM (supersymmetric Yang- Mills). 
Fradkin and Tseytlin considered a quasi-abelian covariantly constant SU{2) 
background as in (3.3) in the 4 spacetime dimensions, with Aj = in the 
remaining 6 dimensions. Then the one loop effective action has the form 

where the important observation is that the form of the integrand ^{s) is 
independent of the dimension d, with the dimension dependence entering in 
the measure of the propertime s integration. Thus, one can compute the 
d = 4 SYM case beginning from d = 10. Analogous to the Heisenberg-Euler 
definitions (1.6), Fradkin and Tseytlin defined the (Euclidean) invariants 



Fl^ = Ji ± ^ J2 - Jl (3.11) 

where Ji = ^F^j^yF^y, and J2 = \F^yF^y. Then for spin 0, i, and 1, one 
finds that for gauge group SU{2): 



sinh(si<'i) sinh(sF2) 
$(Ai/2) = 2['^/2] cosh(sFi) cosh(sF2) $(Ao) , 
$(Ai) = [d + A{smh^{sFi)+smh^{sF2))] $(Ao) . (3.12) 
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where Ao^i/2,1 ^^^ the corresponding operators for fields of spin 0, ^, and 1. 
Thus for d = 10 we obtain 

^s) = -2cI>(Ao) - ^^(Ai/2) + $(Ai) 

= 4(coshsFi - coshsFa)^ $(Ao) . (3.13) 

Therefore, in the dimensionally reduced d = 4 theory one has [106] 

This expression (3.14) is UV finite, owing to the supersymmetry. It is, 
however, IR divergent unless F = F, m. which case -Fi = F2. This IR 
divergence is due to the anomalous magnetic moment term in the spin 1 
operator Ai, as has been noted already in our discussion of QCD effective 
Lagrangians in Section 3.1. 

Another illustrative supersymmetric application of Heisenberg-Euler 
techniques is to the computation of the one loop effective Lagrangian for 
4 dimensional M = A SYM in a constant J\f = 2 SUSY background [107]. 
This M =2 background is characterized by a vector multiplet W = {VFq,, $}, 
and if W satisfies the constancy conditions 

F>aWf3 = D^^Wf^-^ = constant , <& = constant , (3.15) 

then the superfield analogue of Schwinger's proper-time formalism [105] can 
be applied to compute the one loop effective action [107] 

(3,16) 
Here uj{x,y) is the trigonometric function 

/coshx— 1\ /coshy— 1\ / x^ — y^ \ 1 
\ x'^ J \ y'^ J \ cosh X — cosh y y 2 

and ^ and ^ are the fields 

i'^^A.D2(W^\ ^ *., 152 ('«;!), (3.18) 

4$2 V ^ y 4$2 \^ <j2 y V ; 

Notice the clear similarity to the basic Heisenberg-Euler structures in the 
integrand. Supersymmetry combines contributions from different spin fields 
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with particular weights. Buchbinder et al show [107] that this result can be 
made manifestly N = 2 superconformal invariant to yield 

(3.19) 

Here ^ is the full superconformal field, and the function Q{x,y) is related 
to the function uj{x,y) as follows [107]: 

oo 
If Oj{x,y)= Y. Cm,nX^"'y^^ 
m,n=l 

oo 

then n(x,y)= V ^"'" ,, x^ V" • (3-20) 

^^^(2m)(2m + l)(2n)(2n+l) ^ ^ 

The interesting thing about this computation is that the fact that the the- 
ory is superconformal allows one to use the classification of superconformal 
invariants in order to go beyond the Heisenberg-Euler constant field approxi- 
mation. This is an example where added symmetry, combined with the basic 
Heisenberg-Euler structure is a very powerful computational tool [105]. 

3.3. Instanton background 

An important example of nonabelian effective actions concerns the spinor 
and scalar effective action in an instanton background. This was first com- 
puted in the small mass limit by 't Hooft [108]. While not strictly-speaking 
of Heisenberg-Euler form, there are many parallels, since the spectra of the 
massless Dirac and Klein-Gordon operators are known in closed-form for an 
instanton background, as are the associated Green's functions. Nevertheless, 
unlike the Heisenberg-Euler cases, there is no known exact solution for the 
one loop effective action in an instanton background for arbitrary quark (or 
scalar) mass. Consider the case of an SU{2) single instanton background: 

A,{x)=A<^/-^ = ^^h^!^^ (3.21) 

^2 g x'^ + p'^ 

where p denotes the scale of the instanton, r" are the Pauli matrices, and 'q^ya 
are the 't Hooft symbols, which mix the spacetime and group indices [20,108]. 
The corresponding field strength is 

i* pLuyX) = -r-^ ; 7y-^ ■ \6.Z1) 
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Due to the self-duality of the background field, the Dirac operator has a 
quantum mechanical supersymmetry [108-111], which means that the spec- 
trum of the Dirac operator coincides with that of the corresponding Klein- 
Gordon operator, except for a multiplicity factor of 4, and the presence of 
zero modes for the Dirac operator. The net effect is that the the one-loop 
spinor and scalar effective Lagrangians in such a background are related by: 

:« . =-2/:« +U^' --^- ^"^' 



^spinor — "^ ^scalar "^ 2 i 32tt'^ ^^^^ j ^^ ( „2 ) ' (3.23) 

2 

Here u is a renormalization scale mass, and -tSt^ F";,,F";,, is the number den- 
sity of zero modes, with the total number of zero modes being the instanton 
number 

^0 = ^ [d^xF^^F^, , (3.24) 



32vr2 J "■ •- - M^- M^ 

where A^o = 1 for the single instanton field strength in (3.22). The implica- 
tion of (3.23) is that one can compute the (somewhat simpler) scalar effective 
Lagrangian and deduce the corresponding spinor one. 

Using a combination of the solubility of the zero mass Green's functions 
and spectra, and spacetime-dependent mass regulators, 't Hooft computed 
the leading small mass behavior of the scalar effective action as [108] 

silL-lHmp)+a[\)+... (3.25) 

where the constant a( i ) is given by 



«(^) = ^ (7 + In vr) - ^ - i^ « 0.145873 . (3.26) 

Carlitz and Creamer [112] computed the next order term in the small mass 
expansion, and Kwon et al [113] the next-to-next order term, so that we now 
have the small mass expansion 

silL ~ I Hmp) +^{l)+l i^pf iH^p) + 7 - In 2) + O {{m pf) . 

(3.27) 

Gross et al [114] extended the 't Hooft style determinant calculation to 
the important case of finite temperature instantons, or calorons. Very re- 
cently [115], the analogous computation has been performed for the more 
complicated situation of a background caloron field with nontrivial holon- 
omy. 
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The large mass expansion of the one loop effective action in an instanton 
background can be computed using various techniques, ranging from the 
operator product expansion and the derivative expansion [19,20], to the 
heat kernel [113] and inverse mass expansion of [69]. The first two terms 
were computed in [19], and the next two in [113]: 

(1) 11 17 1 232 1 7916 1 

scalar ~ "7^ (mp)2 ~ 735 {m p)^ ^ 2835 {m p)^ ~ 148225 {mp)^ +■■■ ■ 

(3.28) 

A smooth interpolating approximation between the small and large mass lim- 
its (3.27) and (3.28), based on a modified Pade approximation, was proposed 
in [113]. Such an interpolation is of physical interest for the extrapolation of 
lattice results, obtained at unphysically high quark masses, to lower physi- 
cal quark masses. Details of this type of interpolation were studied for the 
soluble Heisenberg-Euler case in [116]. 

4. Two-Loop Heisenberg-Euler effective Lagrangian 

All the previous results have been for the one-loop effective Lagrangian, 
which neglects virtual lines inside the single scalar or spinor loop. In this 
Section I summarize what is known about the two-loop effective Lagrangian. 
In principle, the computation of the two-loop Heisenberg-Euler effective La- 
grangian in QED is completely straightforward, as we only need to compute 
a single vacuum diagram (see Figure 3) with an internal photon line and a 
single fermion (or scalar) loop, where these spinor or scalar propagators are 
in the presence of the background field. These background field propagators 
have been known in closed-form for a long time [12, 117, 118]. 




Figure 3. The two loop diagram for the two loop effective Lagrangian. The double line refers to 
a propagator in the presence of the constant background field, while the wavy line represents the 
internal virtual photon. 



The problem is, however, not completely straightforward because at two- 
loop we need to perform mass renormalization in addition to charge renor- 
malization. Nevertheless, in seminal work in the mid-1970's, Ritus found 
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exact integral representations for the fully renornialized two loop Heisenberg- 
Euler effective Lagrangian in both spinor [88] and scalar [89] QED. These 
were impressive computations, but unfortunately the answers are rather 
complicated double-parameter integrals [see (4.9) and (4.26) below]. I first 
review the structure of Ritus's solutions and then discuss some recent sim- 
plifications that have been found. 

4.1. Two loop spinor QED Heisenberg-Euler effective 
Lagrangian 

Using Feynman gauge for the internal photon line in Figure 3, the two-loop 
spinor QED effective Lagrangian is [88,90] 

2 ,■ 

>cgnor = i y y d'^' tr bM^, x')j^G{x',x)] D{x - x') 

= 4lwr'l Wf'' [>G(P)7^G(,)] D{p - q) (4.1) 
where the position and momentum space propagators are related by 

G{x, x') = e^*(^'^') / -^ e^P-(^-^') G{p) . (4.2) 

Here $(x, x') is a gauge dependent phase factor coming from the line integral 
of the gauge field from x to x'. For a constant background field strength -F^jy, 
using the Fock-Schwinger gauge (2.43), one finds [88,90] 



l2 



m - '-iM" + eF^z, 



tj^ spinor (,2^) X ) ^ —I — -^ / 

exp I -im\s - L{s) + '-z^p^'''{s)z, + '-^a^,F^"'s\ (4.3) 

where we have defined z^ = x^ — x'^, and /3^;y and L are functions of the 
proper time s involving trigonometric functions of the 4x4 constant field 
strength matrix F^j^ as follows: 

P,,{s) ^ [eF coth{eFs)l^ , L{s) ^ hv In (^^^^^^) • (4.4) 

Note that the contour of the proper time s integration is defined to pass 
below the singularities on the real axis, in accordance with the ie prescription 
w? ^ m?' — ie [90]. 

It is now straightforward, but somewhat tedious, to do the Dirac traces 
and spacetime integrals in (4.1) to obtain the following expression for the 



46 Gerald V. Dunne 

bare two-loop effective Lagrangian as an integral over the proper time pa- 
rameters s and s' of the two spinor propagators: 



spinor 1287r^ Jq Jq sin(erys) sin(erys') sinh(ees) sinh(ees') 

xUm^{SS' + PP')lo-il\. (4.5) 

Here (following Ritus's conventions [88]), the two eigenvalues of F^^, are 
called 7] and e, rather than b and a respectively, as used by Heisenberg and 
Euler (1.5). Thus, in a frame in which the magnetic and electric fields are 
parallel we can take rj = B, and e = E. The functions 5, S", P, P', I and 
/q appearing in the two-loop bare effective Lagrangian (4.5) are defined as: 



S{s) = cos(erys) cosh(ees) 
S' = S{s') 

— a a 

a = ei] (coi{ei]s) + cot(erys')) 

2(er/)^ cosh(ee(s — s')) 



sin(e?7s) sin(er7s') 



P{s) = sin(er/s) sinh(ees) 
P' s P(s') . 



/=(i^i„^-a 



q _ p 

_ _ a 

2 



aY a \b — a 
ee (coth(ees) + coth(ees')) 
2(ee)^ cos(er/(s — s')) 
sinh(ees) sinh(ees') 



(4.6) 



The bare effective Lagrangian (4.5) has several different types of divergences. 
These must be regulated and then interpreted physically. Proper time regu- 
larization involves cutting off each proper time integral at a lower bound sq, 
as was done by Ritus [88]. The divergences can then be isolated systemati- 
cally. First there is a subtraction of the field free case, so that £(^' vanishes 
when the background fields vanish. Then there is charge and wavefunc- 
tion renormalization, just as for the one-loop effective Lagrangian, which 
involves identifying a divergent term in C^"^' of the form of the zero-loop 
Maxwell Lagrangian. This is done simply by expanding the integrand to 
quadratic order in the fields rj and e. This divergence can be absorbed by 
redefining the electric charge and the fields as 

ryl/2 -7-1/2 -7-1/2 . . ry. 

e = eoZ^' , r] = riQZ^ ' , £ = £0^3' (4.7) 

where Z3 is the usual multiplicative renormalization factor [119]. With these 
redefinitions, we can express C^"^' in terms of the renormalized charges and 
fields. In fact, note that the combinations er/ and ee are renormalization 
invariant. Third, and more complicated, is mass renormalization. Even 
after dealing with the charge renormalization divergence, there remains a 
logarithmic divergence associated with taking one or other of s and s' to 
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0, while keeping the other one fixed and nonzero. This divergence can be 
separated and one finds that it has precisely the correct form for the one-loop 
renormalization of the electron mass: 



mji 



4^(-l) 



rrir, + 6m 



C 



(1) 



mn 



+ bm 



dC 



R ' 



TTln 



drriQ 



(4.8) 



Note that the second term in (4.8) is of order a^, because 6m'^ and C^ are 
each of order a. Thus it is natural to see this term arising in the compu- 
tation of the two loop effective Lagrangian. The identification of the mass 
renormalization is quite an involved manipulation, and we refer the reader 
to [88,90] for details. In particular, finding the correct finite part of the mass 
shift (5m^ for the relevant renormalization scheme is not straightforward. It 
relies on an independent computation of 6m?, or can be done a posteriori 
by studying the imaginary part of the two-loop effective Lagrangian in an 
electric background, as is discussed in more detail below in Section 4.1.1. 
The final answer for the renormalized two loop effective Lagrangian is 



C 



(2) 

R, spinor 



le 
64^ 

64^ 



ds I ds < K{s, s 



K, 



olsj 



(isi^o(s) { \n{im s) -|- 7 

6 



(4.9) 



where 7 ~ 0.577... is Euler's constant (note that Ritus calls Euler's constant 
In 7), and the functions K{s,s') and Kq(s) are 



K{s,s') 
1 



^-im^s+s') f (M 



2^ee)2 



ss'{s + s') 



Am 



2i 



s + s' 



ppi 

+ 



[Am^{SS' + PP')lQ-il] 



e"'(ry^ 



{2m {ss - 
biss' 



2s^ - 2(s' 



/\2\ 



s + s' 



Kois) 



Am + i 



d_ 

ds 



(er?)(ee) 



tan(er/s) tanh(ees) s 



1 e^ir? 



^2^^ 



(4.10) 



Ritus's expression (4.9) for the renormalized two loop Heisenberg-Euler ef- 
fective Lagrangian is finite and expressed in terms of renormalized quantities. 
It is, nevertheless, a very complicated expression, and it is significantly more 
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difficult to extract physical information from it, along the lines of the appli- 
cations discussed in Section 1.2. For example, one can extract the two-loop 
light-light scattering low energy effective Lagrangian by expanding Ritus's 
answer to quartic order in the fields. One finds [88] 



spinor 



647r4 



m^ 



81^' 



.V + l(,.)^ 



(4.11) 



which should be compared to the one loop result (1.9). Similarly, in the 
limit of a strong magnetic field, one finds [88] 



^P'"""" 1287r4 



In 



er] 

Tim? 



+ constant 



+ 



(4.12) 



from which we can extract the two loop spinor QED /3-function, as discussed 
below in Section 5. This leading strong-field behavior (4.12) should be com- 
pared with (1.54). However, the general weak and strong field expansions of 
Ritus's result (4.9) are very cumbersome, due to the double integral structure 
and the complicated integrand. In Section 6 we shall see that for a self-dual 
background the corresponding expression simplifies dramatically and much 
more can be said about the weak and strong field behavior of £'^'. 



4.1.1. Two loop spinor case for magnetic background 

If we specialize Ritus's result (4.9) to a constant magnetic background, by 
taking rj = B and e = 0, then C^'^' simplifies a bit, but not much. Then one 
finds (after rotating the contours and changing integration variables) 



C 



(2) 

R, spinor 



e^m^ f eB 



(47r)*^ \m? J Jq s' 



OjS _^2 ^ 

e 



'"^/(^■^) / du [L{s,u)-2s'^ 



+ 



u(l - u) V sinh^ 



-|- s coths 
s I 



\ 



-12 



e m eB f^ d^ rr?s/{eB) 



7 — log 



(47r)4 TV? Jq 



1 s 

coths 

s 3 



rn^s 



eB 



+ 



eB 

m?s 



(4.13) 
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where 



L{s,u) = scoths 



^°S ( G{u,s) 



Fi + 



Fo 



[u(l-n)-G(n,s)]2 ' G{u,s)[u{l - u) - G{u,s)] 

F3 



+ - 



G{u, s 



u{l — u)[u{l — u) — G{u, s)] 
coshs — cosh((l — 2u)s) 



2s sinhs 
Fi = 4s(coths — tanhs)G(n, s) — 4ti(l — u) , 

F2 = 2(1 - 2u f^^ ^^ ~ ^'^' + s(8tanhs - 4coths)G(n, s) - 2 , 
smns 

F3 = 4u{l - n) - 2(1 - 2^)!1!M1^^!M _ 4s G{u, s)tanhs + 2 . 

smns 



(4.14) 



While simpler than the general expression (4.9), this magnetic background 
two loop effective Lagrangian (4.13) is still not particularly simple. By a 
direct expansion of the integrand one can compute a number of terms in the 
weak-field expansion [62,120]: 



C. 



2 4 T 

(2) e^m^ 1 

spinor ~ (4^)4^ 



641 4 



eB\ 



^ 135308 /eS^^ 
+ 



1219 



25 \m^J 1225 V"^^ 



791384 /e5\^° 8519287552 /e^^ ^^ 
+ 



1575 V m'^ 



2223375 \ m'^ 



(4.15) 



The coefficients of this expansion alternate in sign and grow factorially fast, 
and the Borel summation properties of the series have been studied in [120]. 
At present, no closed formula is known for these expansion coefficients. 

For an electric background, /^sp-noj. is given by the replacement B^ ^ 
—E"^. Lebedev and Ritus [121] analyzed the singularity structure of the 
integrand in (4.13) to deduce that including the the two loop contribution 
led to: 






-sp 



-sp 



87r3 



fc=i 




(4.16) 
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where a 



A-K 



is the fine structure constant, and the coefficients c^ are 



ci = 



Cfc 



fe-i 

E 



2VkJTi^/l{k-l) 



k>2 



(4.17) 



Note that the two loop contribution has the same exponential form as at one 
loop, but with an overall factor of an, and with a (perturbative) expansion 
in the electric field strength in the prefactor. Ritus and Lebedev were able to 
extract information about the coefficient Ck of the first term in this prefactor. 
For k > 2 the prefactor expansions in (4.16) begin with terms that are 
singular in the limit of vanishing field E ^ 0, which seems potentially un- 
physical. In [121] a physically intuitive solution was offered to this apparent 
dilemma. Their proposal is that if one could take into account all contribu- 
tions from even higher loop orders to the prefactor of the k-th exponential, 
then one would find them to exponentiate in the following way. 



fc2 



+ aTT Ku 



^^U... 



m- 






exp 



kmn 
eE 



fc2 



exp 



kTTml{k,E) 



eE 



. (4.18) 



Thus, it should be possible to absorb their effect into a field-dependent shift 
of the electron mass. Using just the lowest order coefficients (4.17) in the 
small E expansion of Kj^{^), this mass shift reads 



m^{k,E) 



m-\ — akck' 
2 



r-^ I , eE 

' eE ak — 

2 m 



(4.19) 



As shown in [121, 122] these contributions to the mass shift have a simple 
meaning in the coherent tunneling picture : the negative term can be inter- 
preted as the total Coulomb energy of attraction between opposite charges 
in a coherent group; the positive one, which is present only in the case 
k > 2, represents the energy of repulsion between like charges. The pic- 
ture is consistent with the work of Affleck et al [123] who showed that in 
the weak field limit, the dominant result of all higher loop corrections to 
the leading one-loop pair production rate (1.10) is a multiplicative factor of 

It is important to note that this interpretation of the mass shift requires 
the mass m on the right hand side of (4.19) to be the physical renormalized 
electron mass of the vacuum theory. Only in this case does the expansion of 
Kk{^) have the form indicated in eqs. (4.16) and (4.17). It is an interesting 
corollary of the Lebedev-Ritus analysis that the physical electron mass can 
be recognized from an inspection of the two-loop effective Lagrangian alone, 
without ever considering the one-loop electron mass operator. Thus, one 



Heisenberg-Euler Effective Lagrangians 



51 



can fix even the finite part of tlie mass renormalization by demanding tliat 
the physical mass be that which appears in the leading exponential of the 
imaginary part of the effective Lagrangian in an electric field [121]. 

Later, in a Borel analysis [120] of the large order perturbation theory be- 
havior of the two loop result (4.13), a numerical value for the next coefficient 
was extracted, but only for the leading single-instanton term: 



Imf£g)+£g) 



87r3 



1 + avr 1 - 0.44 




+ ... 



eE 



(4.20) 



This numerical approach was not sensitive to higher order terms in the in- 
stanton expansion [the sum over k in (4.16)], as such terms are exponentially 
suppressed. 

Dittrich [124], and Dittrich and Renter [5] also computed this two-loop 
effective Lagrangian, for the case of a magnetic background, and their final 
answer was in a slightly different form from Ritus's result. The strong field 
limits agreed, but it was later found that the weak field expansions did not 
agree, and the cause was found to be a difference in the finite part of the 
mass renormalization [125]. This latter approach [125] used the worldline 
formalism and dimensional regularization, which was also applied to the case 
of a general constant electromagnetic background in both spinor and scalar 
QED in [126]. In Section 7.3 a much simpler, almost trivial, method of mass 
renormalization is presented, for the case of scalar QED. 



4.2. Two loop scalar QED Heisenberg-Euler effective 
Lagrangian 

Ritus also computed [89] the two loop renormalized Heisenberg-Euler ef- 
fective Lagrangian for scalar QED, generalizing Weisskopf's one loop result 
(1.8). For scalar QED, the scalar propagator in a constant background field 
-F^jy is simpler in form than the corresponding spinor one (4.3), but the 
vertices are more complicated as the momentum dependence of the vertex 



changes from p^ to (p^ 
in coordinate space is 



eA^) in a background field. The scalar propagator 



-i%XaF^'•'x' 



^scalar 1,2^) X 



(4vr) 



^ ds 
-7T exp 



-im s — L{s) + -zf3{s)z 



(4.21) 



where (3^v{s) and L{s) are the same functions defined in (4.4) for the spinor 
propagator. Using Feynman gauge for the internal photon propagator, the 
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two loop effective action is [89] 

€L -^^^''^/^'Vc,. - .') {(.|n,.G|.'>(.1n„G|.> 

+ {x\U^GU^\x'){x\G\x') + Ai6{x - x'){x\G\x')} (4.22) 

where we have dropped a term (x|n^G|a;)(a;'|n^G|x'), which vanishes. The 
delta function term in (4.22) corresponds to the tadpole diagram, a new 
feature of scalar QED compared to spinor QED. Using the basic matrix 
elements [89] 



{x\U^G\x') 

/'^ els 1 ( i ^ 



^-i^x^F^-x'^ /.oo ^g ^ 



(47r 

{x\GIi^\x') = 
, e-^f ^M^^'"'^- r^ ds 1 



, 29 [^^'^ ~ ^^^''l ^'^ ^^P > ~^^ ^ ~ ^^''^^ ^ 7^/^(^)^ 



(x|n^Gn^|x') = -i5{x - x) - m?{x\G\x) + —{zFFz){x\G\x') (4.23) 



=,2 

2" 

one finds that the bare, unrenormalized, two loop effective Lagrangian for 
scalar QED is 

,(2) ie^ f°^ f°° , (er/)2(ee)2e^*'"'(^+*') [m^/o + 1/] 



/::,,. = — -^ / ds / ds 



scalar 2567r'^ Jq ^0 sin(er7s) sin(er/s') sinh(ees) sinh(ees') 

(4.24) 

where the forms of the functions /q, /, a and h are the same as in the spinor 
case (4.6), but the functions p and q entering in the definition of / are 
changed from those in (4.6) to 

p = li&Tif' (cot(erys) cot(er/s') + 3) , 

q = 2{eef (coth(ees) coth(ees') - 3) . (4.25) 

The bare effective Lagrangian (4.24) also contains several divergences, which 
may be regulated using a proper time cut-off sq. These divergences corre- 
spond, as in the spinor QED case, to charge and mass renormalization. 
Proceeding as in the spinor case, one finds [89] 

(2) _ ^^' r..r.,ir.^.i. Kois) 



""R, scalar -i 9877-4 
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(4.27) 



Just as in the spinor case (4.9), the general weak and strong field expan- 
sions of Ritus's scalar QED result (4.26) are very cumbersome, due to the 
double integral structure and the complicated integrand. Also, note that it 
is very different from the spinor result, in the sense that the integrands are 
very different. We return to this comment later when we consider self-dual 
backgrounds in Section 6. 

Nevertheless, it is still possible to extract from (4.26) the leading weak- 
field and leading strong-field behavior [89]. For weak fields, expanding the 
integrand to quartic order in the fields. 



scalar 



275 
2592 



(r/^ 



^^)^ + ^(^^) 



+ 



(4.28) 



which should be compared to the one loop result (1.32). Similarly, in the 
limit of a strong magnetic field. 



.(2) 
"scalar 



128^ 



In 



erj 
TTm? 



+ constant 



(4.29) 



from which we can extract the two loop spinor QED /3-function, as discussed 
below in Section 5. Indeed, the fact that (4.12) and (4.29) have the same 
leading coefficient tells us that the two loop /3-function coefficients in spinor 
and scalar QED are equal [see (5.13) and (5.14) below]. 



4.2.1. Two loop scalar case for magnetic background 

The general expression (4.26) can be specialized to just a magnetic back- 
ground by setting e = and r] = B, and one obtains an expression similar in 
form to (4.13) [62]. This expression is, however, still not particularly simple. 
By a direct expansion of the integrand one can compute a number of terms 
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in the weak- field expansion [62]: 
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(4.30) 



The coefficients of this expansion alternate in sign and grow factorially fast, 
just as in the spinor case discussed in the previous section. At present, no 
closed formula is known for these expansion coefficients. 

To conclude this section on the two loop effective Lagrangians, we briefly 
mention that the world-line techniques have been applied to two loop effec- 
tive Lagrangians in Yang- Mills theory [127], where there are infrared diver- 
gences due to the masslessness of the gluon fields. Kuzenko and McArthur 
[128] have computed the two-loop Heisenberg-Euler effective Lagrangian for 
N = 2 SUSY QED, finding a result expressed as a double parameter integral 
similar in structure to those in the spinor and scalar QED results of Ritus. 

5. /3-functions and the strong-field limit of effective 
Lagrangians 

In quantum field theory there is a close connection between the short- 
distance behavior of renormalized Green's functions and the strong-field limit 
of associated quantities calculated using the background field method. This 
correspondence leads, for example, to a direct relation between the perturba- 
tive /3-function and the strong-field asymptotics of the effective Lagrangian, 
as was first established by Coleman and Weinberg [87]. Ritus developed 
this relation for QED using the renormalization group, with the assumption 
that the strong-field limit of the renormalized effective Lagrangian is mass- 
independent [88-90]. Another, equivalent, derivation which invokes the scale 
anomaly [129, 130] in a massless limit, has been given in various forms by 
many authors [8,86,131-133]. This derivation is reviewed below for QED in 
Section 5.1. 

I stress that computing the /3-function, which characterizes the scale de- 
pendence of the running coupling, is a much simpler problem than finding 
the full renormalized effective Lagrangians, (4.9) and (4.26), as one only 
needs the leading strong field behavior of C In this approach, one uses the 
external field as a probing scale, instead of the external momentum of a 
self-energy diagram (which is the usual text-book approach). One adapts 
the Gell-Mann Low [16] renormalization group argument for n^(j(q^,/i^) at 
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large q^, to the effective Lagrangian C{eF, fi'^) to relate the strong-field limit 
of C to the /3- function coefficients. A similar idea was used at two loop by 
Shifman and Vainshtein to compute the QED /3-function using the operator 
product expansion [134], and has recently been extended to the three-loop 
Yang- Mills /3-function [135]. There is an immediate combinatorial advantage 
to this background field approach, as there are many fewer diagrams at a 
given order, since there are no external lines. Furthermore, each diagram 
has fewer vertices and propagators. The price, of course, is that the spinor 
(or scalar) propagators are not free ones, but are in the background field. 

Here I concentrate on QED, but these ideas generalize, and have been 
applied to, nonabelian theories. Begin by recalling some basic facts about 
the QED /3- function [119, 136]. The /3- function encapsulates the scale de- 
pendence of running coupling, and is defined as 

2 

where the natural expansion parameter turns out to be a = ^, with a = |;p 
being the fine structure constant. At present, the spinor QED /3- function is 
known to four loop order. Beyond two loop order the /?- function coefficients 
are scheme dependent; for the MS scheme the four loop result is [137]: 



/?Ms(a) 



a^ + 4a^ - ^a^ 
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a^ + ... 



(5.2) 



It is worth emphasizing how difficult these computations are. The two loop 
result, Aa^, was obtained in 1950 by Jost and Luttinger [138], the three loop 
result by Rosner in 1966 [139], and the four loop result in 1991 by Gor- 
ishny et al [137]. (One is tempted to extrapolate that the five loop result 
will require 36 years to compute!) The four loop result was computed us- 
ing computer programs to organize and manipulate the many hundreds of 
Feynman diagrams involved. This complexity contributed to the earlier pub- 
lication of a different result for the four loop term by the same group [140], 
this result being retracted four years later after errors were found in some 
of the computer subroutines. This illustrates how difficult these calculations 
are, and motivates other ways to understand and perform the computations. 
Nevertheless, the techniques [137] used to produce this landmark result (5.2) 
for the four loop QED /3-function coefficient have now been applied to QCD 
where we now know the /3-function and the anomalous quark mass dimension 
to four loop order [141], a truly remarkable achievement. 

Another /3-function of interest is the so-called " quenched" /3-function (or 
Fi function), in which one only includes a single fermion (or scalar) loop, 
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neglecting self energy corrections to internal photon lines. This /3-function 
is scheme independent, essentially because the corresponding single-fermion- 
loop diagrams only diverge like a single logarithm. At four loop order, the 
quenched spinor QED /3-function is 

pQ{a) = -a^ + Aa^ - la"^ - A6a^ + . . . (5.3) 

o 

which has the remarkable property that the coefficients are rational num- 
bers. This is despite the fact that at intermediate stages of the computation 
individual diagrams, and groups of diagrams, contribute irrational terms 
such as C(3) and C(5), but these all cancel at the end when the quenched 
/3-function is computed. This fascinating property was already noticed at 
the 3-loop level by Rosner in 1966, and has been studied in detail over the 
years [142]. Despite impressive recent progress by Broadhurst et al [143] who 
have identified the cancellation of such irrationals with a relation between 
Feynman diagrams and knots, a truly deep understanding of this behavior 
is still lacking. A promising approach is the worldline formalism [62] which 
combines the relevant Feynman diagrams into a much more compact repre- 
sentation, making the large cancellations somewhat more 'natural'. This has 
been studied at the two-loop level [60], but little is known beyond two loop 
at present. Even at a more basic level, we do not know how the perturbative 
series for the /3-function behaves. It is presumably divergent, but we do not 
know if it is Borel summable. The first four signs of — , — , -|-, -|-, do not give 
us much basis for predicting the alternating or nonalternating behavior of 
the full series. These comments motivate the exploration of various different 
ways to compute /3- functions, and the strong field limit of a Heisenberg-Euler 
effective Lagrangian technique gives a particularly efficient way to do such 
a computation. 

5.1. General Argument 

In this section I review the general argument [8,86,88,90,131-133,144] relat- 
ing the strong-field asymptotic behavior of the effective Lagrangian to the 
perturbative f3 function. I present the argument for QED, but it is more 
general. Consider an abelian gauge field coupled to spinor or scalar matter 
fields, which are either explicitly massless or which have a well-defined mass- 
less limit. Then the trace anomaly for the energy-momentum tensor states 
that [129, 130] 

W = §r|J(^M^)^ (5-4) 
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where e is the running coupHng, and /3(e) is the /? function, defined below 
in (5.8). The expectation value of the energy-momentum tensor can also be 
related to the effective Lagrangian for a constant background field strength 
F ■ 

{en = -v'"'C,s + 2^ . (5.5) 

These two relations, (5.4) and (5.5), determine the effective Lagrangian to 
be of the form 



1^ 

4 e2(t) 



'Ceff = -7i277T^M^^'^^ (5-6) 



where the " renormalization group time", t, is expressed in terms of the scale 
set by the field strength, serving as the renormalization scale parameter 

IJ? ~ e\F\, 

1 , /e2|F2|- 



and ^0 denotes a fixed reference scale at which, for example, the value of the 
coupling may be measured. Note that in this argument the field strength 
plays the role which is usually associated with a momentum transfer Q^. 
This already suggests at a very basic level why the strong-field and short- 
distance limits are related. 

The P function is defined in terms of the running of the coupling as 

mt)) = ^ ■ (5.8) 

To see how this solution (5.6) leads to an explicit connection between the 
strong-field asymptotics of £cff and the perturbative /? function, note that 
(5.8) can also be expressed as 

^(*) de' 

where e = e(0). Making a perturbative expansion of the [3 function 

/3(e) = /?ie3 + I32e' + ... (5.10) 

the relation (5.9) determines the running coupling, e(t), in terms of e as 

^-- = -^-2l3it-2P2e^t + 0{eH'^) . (5.11) 
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Inserting this into (5.6), the strong-field asymptotics of the effective La- 
grangian is, to two-loop order, 

£eff ~ ^ (2/3ie2 + 2/326^ + • . . ) F^.F^^'^ In f ^^) , (5.12) 

where, as is conventional, we have subtracted the classical Lagrangian, 
— \F^i,F^^ , from >Ceff. Note the appearance of the /3-function coefficients 
/?! and (32 in the prefactor of the strong- field behavior of the C. At higher 
loops the structure is slightly more complicated, but knowledge of the leading 
strong-field behavior up to a given loop order determines the /3-function coef- 
ficients to that order, within the corresponding renormalization scheme [90]. 
In order to illustrate this correspondence explicitly, we recall that the 
QED fj functions, for spinor and scalar QED, to two-loop order, are 



e^ e' 



/3spinor- i27r2 + 647r4 + ■ • • ' ^^-^^^ 



o3 ^5 



/3.cala. = ^ + ^ + .... (5.14) 



5.2. Explicit example: constant magnetic field background 

Equation (5.12) gives a direct correspondence between the perturbative (3 
function coefficients and the strong-field behavior of the effective Lagrangian. 
We now compare this with some explicit results where the effective La- 
grangian is known. First, consider the Euler-Heisenberg effective Lagrangian 
for a constant background magnetic field, of strength B. 

At one loop, the leading strong-field asymptotics is [recall (1.54) and 
(1.63)]: 

^(1) magnetic e B , / ^S \ /K 1 r:^ 

4pi„o, ~^j?'"(;;?j + '-- . (5-15) 

4r-"=~Sin(:i)+.... (5.10) 

Noting that —^F^^F^^ = —^B"^, and comparing with the correspondence 
(5.12), we deduce that /3^''"^°'^ = j^ and f3f^^^^^ = j^, in agreement with 
the one-loop (3 function coefficients quoted in (5.13) and (5.14). 

At two loop, the corresponding leading behaviors are deduced from Ri- 
tus's results. Recalling the leading strong field behaviors (4.12) and (4.29) 



we have 
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C^'') '^^s°^"^ ^ ^-^ In (— ] + (5 17) 



,(2) magnetic 



,(2) magnetic & B , f eB 



£;ei,r-"~I5i^ln(^^j+.... (5.18) 

Once again, comparing with the correspondence (5.12), we deduce that 
^spmor _ ^scalar _ ——^^ in agreement with the two- loop (3 function coef- 
ficients quoted in (5.13) and (5.14). 



6. Two loop Effective Lagrangian in a self-dual background 

Ritus's results for the two loop Heisenberg-Euler effective Lagrangian are sig- 
nificantly more complicated in form than the one-loop results. This restricts 
somewhat the investigation of two loop effects. For example, as discussed in 
Section 4.1.1, we have only very limited information about the two loop gen- 
eralization of the pair production formulae (1.25) and (1.37). In this section 
I review some recent developments in which the two loop Heisenberg-Euler 
effective Lagrangian has been shown to simplify dramatically when the back- 
ground field is self-dual. 

6.1. Closed-form effective Lagrangians in a self-dual 
background 

Consider restricting the constant electromagnetic background to be self-dual 

F^u = F^v (6.1) 

where F^j_y = ^e^^p^FP'^ is the standard dual electromagnetic field strength. 
In QED, which is an abelian theory without instantons, such a classical 
background field is unphysical in Minkowski space as it means E = iB 
(or B = —iE). However, Duff and Isham showed that such nonhermitian 
classical fields arise as vacuum to coherent state matrix elements of hermitian 
quantum field operators [97], which means that effective actions in such 
backgrounds are physically meaningful generating functions. We will discuss 
this further below in terms of helicity amplitudes. Also, Euclidean classical 
configurations can be used as probes of real QED properties, such as the 
/3-function and the anomalous mass dimension. 

When the constant background is restricted to be self-dual, the fully 
renormalized two-loop Heisenberg-Euler effective Lagrangian takes a remark- 
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ably simple closed- form, for both spinor and scalar QED [145-147]: 



C. 



(2) 



spinor 



,(2) 
"'scalar 



vi^[3e^(.)-e'( 



Q 



47^2 



■e{^)-^'{^) 



(6.2) 
(6.3) 



Here / is the field strength parameter, \F^uF^^'^ = /^, and k is the natural 
dimensionless parameter 



m 



K = 



2e/ 



(6.4) 



The ubiquitous function S,{k) is essentially the Euler digamma function 
V(K) = Alnr(«): 



^{k) = —K ( iP{k) — ln(K) + 



2k 



(6.5) 



The subtraction of the first two terms of the asymptotic expansion [22, 23] 
of tp{K) correspond to renormalization subtractions, as will be shown below 
in Section 7. 

It comes as something of a surprise that in this self-dual background both 
parameter integrals in each of the expressions (4.9) and (4.26) for the two 
loop effective Lagrangian can be done in closed form to produce the simple 
answers (6.2) and (6.3). If we simply take Ritus's expressions (4.9) and 
(4.26) and specialize to the self-dual case by taking ij = ie, then it is not 
at all obvious that the resulting integrals can be reduced to such a simple 
closed form as in (6.2) and (6.3). 

It is also interesting to note that in such a self-dual background, the one- 
loop Heisenberg-Euler effective Lagrangians (1.2) and (1.8) for spinor and 
scalar QED are also naturally expressed in terms of this same function ^{k): 



c 



(1) 

spinor 



£ 



(1) 

scalar 



87r2 
87r2 

(47r)2 
(47r)2 



dt 

-IT e 



'2Kt 



2*2 
t^ coth^ t-1- —- 



lnK + C'( 

12 ^ ^ 



-1)+E(k) 



(6.6) 



dt 
t^ 



-2Kt 



sinh t 3 



12 



lnK + C'(-l)+H(K) 



(6.7) 
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Here we have defined the function (which is closely related to the Barnes G 
function [33,34]) 

H(x) = —X In r(a;) -\ In a; \- I dy In T{y) (6-8) 

2 4 2 Jo 

which has the property that H'(k) = S,{k). 

The dramatic simplicity of the two loop results (6.2) and (6.3) raises three 
obvious questions: 

• why are these expressions so simple? 

• why are these expressions so similar? 

• why is the particular function ^{k) so special? 

The answers lie in the three-way relationship between self-duality, helicity 
and (quantum mechanical) super symmetry. 



6.1.1. Self-duality and helicity 

Self-dual fields have definite helicity [96,97]. Indeed, the self-duality condi- 
tion (6.1) is just another way of writing the helicity projection: 



For anti-self-dual fields the other helicity projection vanishes, so that the 
photon field has the opposite helicity. 

It is well-known that scattering amplitudes for external field lines with 
like helicities are particularly simple [148-150]. Since the effective action 
for a self-dual field is the generating function for like-helicity amplitudes, 
it is consistent that the effective action in a self-dual background should 
be simple. However, almost all of these helicity amplitude results are for 
massless particles on internal lines (but see for example [151]), while here 
we see a generalization to massive particles. For example, in massless QCD, 
the tree level amplitudes with all external helicities alike, or all but one alike 
vanish identically [149]: 

M\pi+;p2+;...;pn+] = , M[pi-;p2+; • • • ;Pn+] = . (6.10) 

For massless QED at one loop, amplitudes with all external helicities alike 
vanish except for the case of four external lines, and these are simple [152] 

M[pi+;p2+;...;pn+] = , n/4 

, ie^ (12)* (34)* 
M[pi+;p,+;p3+;P4+] = ^^Y^ • (6.11) 
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Here the RHS is written in terms of a convenient spinor-helicity basis [148- 
150] for the momenta: {ij) = U-{pi)u^{pj), with u±(p) being fundamental 
spinors of definite hehcity. Similarly, when all but one external helicity is 
alike, the amplitudes vanish except for four external lines, where it also has 
a simple expression 

M\pi-;p2+;...;Pn+]=0 , n/4 

ie'^ (12)(34)*(24)* 
M[p,-;p.+;p,,+;p.+l = ^i^jA_iA^ ^ (6.12) 

These massless helicity amplitudes, with 4 external legs, have been extended 
to two-loop level [153], where they also exhibit remarkable simplicities. 

The results (6.2) and (6.3) provide information about amplitudes with 
any number of external lines, with all helicities being equal. Massive QED 
turns out to be rather different from massless QED, but still shows dramatic 
simplicity when considered in a helicity basis. For example, a new prediction 
[146] is that for massive QED the low energy limit of the one-loop and two- 
loop amplitudes with all external helicities -|- behave as: 

rW[fci,e+;/c2,e+;...;A;^,e+] = -2— -— |^2lv^c!v/2 Xn (6.13) 

r(2)[/ci,e|;A;2,e^;...;A;Af,e^] = -Ian ^^,2^2N-4 ^N/2 ^^ ■ (6.14) 

where the momentum dependence of the low energy limit is naturally ex- 
pressed in terms of symmetrized helicity basis elements xn'- 

^^ ^ ^^ {[12]2[34]2 ... [(AT - l)iV]2 + an perms} (6.15) 

where [ij] = {k^\kj) denotes a basic spinor product [148-150]. The numer- 
ical coefficients c„ and c„ appearing in (6.13) and (6.14) come from the 
weak field expansion of the one-loop and two-loop effective Lagrangians for 
a self-dual background [146] : 

(1) _ ^2n 

^" ~~2n(2n-2) 

In massive QED, the [+ -|- -|- • • • -|-] helicity amplitudes for A^ > 6 do not 
vanish at one or two loop, in contrast to the one loop case (6.11) for massless 
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QED. On the other hand, the amphtudes with all but one helicity alike, 
the [— + +•••+] amplitudes, vanish in the low-energy limit in massive 
QED, a result that is true to all orders [146]. These results have since been 
generalized to all helicity combinations (up to 10 external legs) by making 
an explicit helicity expansion of Ritus's general result (4.9). In fact, in the 
low energy limit all "odd -" amplitudes vanish, to all orders [154]. 

6.1.2. Self-duality and quantum mechanical supersymmetry 

Another consequence of the self-duality of the background is that the corre- 
sponding Dirac operator has a quantum mechanical supersymmetry. That 
is, apart from zero modes, the Dirac operator has the same spectrum (but 
with a multiplicity of 4) as the corresponding scalar Klein-Gordon opera- 
tor [108,109,111]. At one-loop this implies [compare to (3.23)] 

4L = -24;L + ^(i{)'>n(^) (0.17) 

where ^o = [^) is the zero mode number density. The logarithmic term 
in (6.17) corresponds to the zero mode contribution, as can be verified di- 
rectly from the one loop Heisenberg-Euler results restricted to a self-dual 
background [144]. Renormalizing on-shell (i.e., /U^ = m^), we find that the 
spinor and scalar effective Lagrangians (6.6) and (6.7) are proportional to 
one another for a self-dual background, in such a way that the supersym- 
metric combination vanishes: 

-^spinor + 2-^ scalar = ^ (6.18) 

These properties are familiar from one loop instanton physics [96, 108, 109]. 
Now consider the implications of self-duality of the background at the two 
loop level. The two-loop the effective action is not simply a log determinant, 
so the situation is more complicated and we do not expect a SUSY relation 
like (6.18) to hold. Nevertheless, the quantum mechanical SUSY of the Dirac 
operator relates the spinor propagator to the scalar propagator via simple 
helicity projections [155], 

5 = -H»-,„)G(i±^VGp(i^) + " + ^*'^'""^ 



m 

(6.19) 

This has the consequence that simply doing the Dirac traces in the spinor 
two loop effective Lagrangian (4.1), one finds that it can be written as the 
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sum of two terms involving matrix elements of the scalar propagator, and 
moreover these are the same two matrix elements of the scalar propaga- 
tor that appear in the scalar QED effective Lagrangian, but with different 
numerical coefficients [144] : 

£(2) = ^J Jj'^ [A{x\D^G\x'){x'\D^G\x)+B{x\G\x'){x'\D^GD^\x)] , 

spinor QED : A = -4 , B = 8 ; scalar QED : A = -1 , B = -1 (6.20) 

This structure explains why the two loop answers (6.2) and (6.3) for spinor 
and scalar QED have such a similar form, involving just two terms with 
different numerical coefficients. 

6.1.3. Propagators in self-dual background 

A more prosaic reason for the simplicity of the two-loop expressions (6.2) and 
(6.3) is that for a self-dual field the square of the matrix F^^ is proportional 
to the identity matrix: 

F^,F,p = -f6^p . (6.21) 

This dramatically simplifies the propagators of spinors or scalars in the back- 
ground field. For example, for a scalar particle, the Lorentz matrix structure 
in (4.21) disappears, leaving (neglecting the unimportant phase factor) 

V47ry Jq smh^(e/t) 

Gscalar(p) = H J' ,, g^^^t-fj ta„h(e/.) _ ^g 33) 

Jo cosh [ejt) 

Note that these propagators are functions of {x — x')'^ and p^ respectively, 
which greatly simplifies computations involving the propagators. Also note 
the duality between the momentum and coordinate space propagators ~ the 
coordinate space form is well suited to a weak-field expansion, while the 
momentum space form is well suited to a strong-field expansion. The propa- 
gators satisfy simple differential equations 

+m)G(p) = l+^-J ^-^ 
{-dl + m^)G{x, x') = 5{x - x') - (^\ (x - x'fG{x, x') . (6.25) 



(p^ + m^)G(p) = l+^^J ^-^, (6.24) 
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The simplification of tlie self-dual propagators is even more dramatic for 
massless particles, where for scalars 

Gsca\iir{x,X ) = - — tt-. -tit , Gscalar(p) = 5 • (6.26) 

47r"'(x — x')^ p"^ 



6.2. Weak field expansions in self-dual case 

The weak field expansions are trivial to derive because we have simple closed- 
form expressions for the two loop effective Lagrangians in terms of the func- 
tion ^(k) in (6.5), whose analytic properties are completely understood, since 
it is essentially the digamma function. To derive the weak field expansion, 

2 

first recall the definition of k as «; = ^r, so that the weak field limit is the 
large n limit. This weak-field expansion follows directly from (6.5), together 
with the asymptotic (large x) expansion of the digamma function [23]: 

^(-)~l---^-f;^ (6-27) 

k=l 

where B2k are the Bernoulli numbers. Then the weak-field (large k) expan- 
sion of (6.2) is 



where the two- loop expansion coefficients are (for n > 2): 



Cinor(-) - -2«-^^ ^^4^^™^ (6.28) 



~ (27r)2 \ 2n - 2 '^'""^ + -^ 2^ 2k {2n - 2k) j ' ^^^^^ 



The leading term in the weak-field expansion (6.28) is 



^(2) ^.^.. e'/' 5 



<pinor('^) - -Zlh^ (6-30) 



agreeing with Ritus's result (4.11) with the self-dual replacement r] = ie = f. 
Similarly, for scalar QED the results are almost the same, but with dif- 
ferent coefficients: 



4 °° 1 

"^ V^ ^(2) scalar -*- 
-scalar V-; '^ "" J^ 2^ C„ -^ 



4cL(-)--7^y:c(^^-^'-4r (6.31) 



n=2 
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where the two- loop expansion coefficients are (for n > 2): 



J2) scalar _ 1 f ^^ ~ ^ g , 3 ^^ g2fc g2n-2fc ) ,^ „„>^ 

'" -(2^\2;rr^^^"-^ + 2^^(2n-2fe)/- ^^'^^^ 

The leading term in the weak-field expansion (6.28) is 

agreeing with Ritus's result (4.28) with the self-dual replacement r] = ie = f. 
Another advantage of having a closed-form expression for the two loop 
effective Lagrangians is that it becomes simple to study the nonperturbative 
contributions to the imaginary part. For example, real k is like the magnetic 
B field case, while imaginary k is like the electric E field case [147]. When k 
is imaginary, using the properties of the digamma function we find the exact 
instanton expansions (here for scalar QED) 



(47r)3 k2 Z^ I Z^ 21k?^ j 



Note that the leading terms in the weak field (large n) limit differ by a 
factor of avr, just as in the electric field case discussed in Section 4.1.1. 
Also, observe that the two loop imaginary part has a prefactor of the same 
form as that proposed by Lebedev and Ritus in (4.16), but in this self-dual 
case we have a closed-form expression for all terms in this prefactor series. 
We see moreover that this prefactor series is an asymptotic series. While 
this self-dual background is unphysical in Minkowski space, these results are 
suggestive of what to expect for a real electric field background at the two 
loop level. 

6.3. Strong field limits in self-dual case 

The strong-field (small n) expansion follows directly from (6.5), together 
with the small argument expansion of the digamma function 



^ oo 

^(x) ^ + y^(_i)fe^(fe)x'=-i . (6.35) 



X 

fc=2 
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This leads to the following strong- field behavior for spinor QED: 
2am'^ 1 



jC^ 



spinor 



(47r)3 K 



3 ^2 



sQ+T^ + z^lnK^ +3ff;(-l)'C(0'^'"M 

-6^+7/^+ Kin/.) (f^i-iyaiW-'j 

oo 

- (1 + 7 + In k) + J^(-l)^/ ClO't'^^ 



1=2 

The leading strong-field behavior for spinor QED is 

^spmor 32^4 \^ ^2 ^ 

Similarly, for scalar QED we find 
,(2) _ am'^ 1 



(6.36) 



(6.37) 



C. 



scalar 



(47r)3 K 



oo 



1=2 



The leading strong-field behavior for scalar QED is 



C 



(2) 



„4 J^2 



t,JM.\ 



scalar 54^4 \^ ^2 ^ 



(6.38) 



(6.39) 



6.3.1. Strong-field limit and (3-function in self-dual case 

Since in some sense the self-dual background is the simplest type of back- 
ground, as it leads to the simplest propagators, this suggests it should be a 
useful probe for computing /3-functions. From the one-loop Heisenberg-Euler 
expressions (6.6) and (6.7), and the closed-form two-loop expressions (6.2, 
6.3) we find 



4:1, ~ig^ 






scalar q^^A 



/^In 



e/ 



m^ 



spmor 



spmor 



^^ f2 ,„ / e/ 



247r2 



/^InRV 



m^ 



e" f2 i„ / e/ 



327r4 



rin 



m^ 



(6.40) 
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Recalling the argument in Section 5.1 relating the coefficient of the leading 
strong-field behavior of the effective Lagrangian to the /9-function coeffi- 
cients, we see that the scalar QED results in (6.40) give the correct one 
and two loop scalar QED /3-function coefficients [compare with (5.12) and 
(5.14), and recall that the Maxwell Lagrangian is /^]. However, for spinor 
QED comparing with (5.13) we see that the correspondence appears to fail. 
The reason is that in a self-dual background spinor QED has zero modes, 
and these preclude the massless limit on which is based the argument in 
Section 5.1 relating the /3-function and the strong field limit of C. This ar- 
gument must be modified to account for the zero modes [144]. At one loop 
the mismatch is due to the infrared divergence of the unrenormalized £spinor, 
as is familiar from instanton physics [108, 109]. But at two loop there is no 
IR divergence in the bare Lagrangian and the zero modes actually enter via 
the mass renormalization [144]. 



7. Diagrammatic approach to two loop effective Lagrangians 

It is natural to ask if the remarkable simplifications of the two loop results for 
a self-dual background might extend to even higher loops. To go beyond two 
loops one should take advantage of the great progress that has been made 
recently in understanding the structure of higher-loop quantum field theory 
(without background fields) . The general strategy is to manipulate diagrams 
to reduce the number to a much smaller set of so-called "master diagrams" 
which need to be computed. This has led, for example, to many new two-loop 
results for QCD scattering amplitudes [156]. Kreimer [157] and collaborators 
have discovered an elegant Hopf algebra structure underlying the seemingly 
messy jumble of ordinary Feynman diagram perturbation theory. I conclude 
this review with some brief comments and speculations about how these 
techniques can be extended to incorporate background fields [158]. 



7.1. Why $ and ^' ? 

The first question to address is why are the simple expressions (6.2) and 
(6.3) for the two loop Heisenberg-Euler effective Lagrangians in a self-dual 
background expressed in terms of the function ^(k) and its derivative, where 
recall that ^ was defined in (6.5) as essentially the Euler digamma function. 
The first hint comes from the following facts that for a self-dual background 
the scalar propagator loops, evaluated using (6.23), turn out to be simply 
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related to the ^{k) function: 




o 



d'^p 



[G{p) - Goip)] 

mp)r - (Go(p)) 



m^ ^(k) 



(47r)2 K 



(47r)2 



(7.1) 
(7.2) 



Here the double lines refer to scalar propagators in the self-dual background 
and the single line is the free scalar propagator, while the dot on a propagator 
refers to the propagator being squared. Thus, ^ and ^ are natural one loop 
traces for the self-dual background. 

Given (7.1) and (7.2), we can write the closed-form expressions (6.2) and 
(6.3) for the two loop effective Lagrangians in diagrammatic form: 



spinor QED : 




-6e' 




+ 



27r2 




(7.3) 



scalar QED 
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47r2 
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(7.4) 



Here the notation is that the triple line loop on the LHS of (7.3) refers to 
a spinor propagator in a self-dual background, while the double- line loops 
[including those on the RHS of (7.3)] refer to a scalar propagator in the self- 
dual background. This shows the remarkable result that the two loop fully 
renormalized answers are expressed naturally in terms of one loop quantities. 
Qualitatively, we can write: 



two loop = (one loop) + (one loop) 



(7.5) 



Interestingly, such a relation with two loop quantities being expressed as 
squares of one loop quantities plus a one loop remainder has been found 
recently [159] in the amplitudes of 4 dimensional super Yang-Mills theory, 
which is a very different theory from QED. This suggests something deeper 
is at work here. This is discussed in the next sections. 
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7.2. Background field "integration-by-parts" rules 

Indeed, we can go further than the quahtative statement (7.5) and derive the 
results (7.3) and (7.4) by algebraic means. First, we identify the source of 
the coefficient factors — 6e^ and |e^ which appear in front of the (one loop) 
terms in (7.3) and (7.4). Note that in free QED {i.e., with no background 
field) it is a straightforward exercise to show that in 4 dimensions 



spinor QED : 
scalar QED : 



-6e' 



-\2 



(7.6) 
(7.7) 



(The loop on the RHS is a scalar loop in each case.) Thus, the two loop free 
vacuum diagram can be expressed in terms of a simpler one loop diagram. 
These results can either be derived by computing each side using dimen- 
sional regularization, or a quicker proof follows from an integration-by-parts 
argument (see below). Notice that the coefficients of the (one loop) parts 
in this free case are exactly the same as the corresponding coefficients in the 
background field case, for both spinor and scalar QED. This is no accident, 
as I now illustrate for the case of scalar QED (for spinor QED the argument 
is similar). 

Consider the derivation of (7.7) using dimensional regularization and 
integration-by-parts [142, 160, 161]. By purely algebraic manipulations it 
follows that 
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(7.8) 



where the dotted line denotes a massless scalar propagator. The first term 
has been written as the square of a one loop diagram but the second term 
is apparently still two loop. However, using simple integration-by-parts ma- 
nipulations, this two loop diagram can also be written as a square of a one 
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loop diagram: 
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Thus, the apparently 2-loop term on the RHS of (7.8) is a square of a one 
loop diagram, leading to 

2 

(7.10) 



©=y 



d-l 



d-3 



This reduces to (7.7) in d = 4, and I stress that this result has been derived 
without doing any integrations, only making simple algebraic manipulations 
on the integrands. 

Now consider the analogous manipulation in a self-dual background [158]. 
First, we extend the scalar propagator in a self-dual background to arbitrary 
dimensions by taking multiple copies of the block diagonal structure of F^y. 
That is, we take as our definition of "self-dual" fields in d 7^ 4 the property 
(6.21), since the notion of self-duality in (6.1) is clearly tied to d = 4. This 
is equivalent to dimensional regularization in the worldline formalism [62]. 
Then the scalar propagator (6.23) becomes 
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dt 
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Then we can repeat the algebraic steps in (7.8) to obtain 
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where we have chosen to isolate this particular coefficient of the square of 
the one loop propagator trace motivated by the free- field result (7.10). 



7.3. Algebraic view of mass renormalization 

The advantage of the manipulation in (7.12) is that it makes the mass renor- 
malization (which was a very difficult part of previous two loop computa- 
tions) almost trivial. To see this we subtract the free field two loop diagram 
from the background field diagram 






+ 0{f 



and then simply complete the square in the first terms: 




n2 
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(7.13) 



(7.14) 



The cross-term in (7.14) is immediately identified with the mass renormal- 
ization because 



6m 
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^d-3 
which can also be derived algebraically. Furthermore, 
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The /^ term in (7.16) contributes to the charge renormalization, and so 
(7.13) can be written as 
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Observe that the first term is now completely finite, so that we can set d = 4, 
and by (7.1) we obtain the first term of the final answer (7.4), the (one loop)^ 
piece, without doing any integrals at all! The second term is manifestly the 
mass renormalization term, and so is absorbed by mass renormalization (4.8). 
The only remaining divergence can be proportional to the bare Maxwell 
Lagrangian f"^, which is then subtracted by charge renormalization. It is 
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simple to isolate and subtract this piece, leaving an 0{f^) term, whose kernel 
in the d ^ 4 limit reduces to a momentum delta function: 

0{f) '=i -A^'e\eff I ^^ [G{p)G{q) - Goip)Go{q)] 5{p - q) 



47r2 



' (ef)' 



e , »^2 




(7.18) 



where in the last step we have used (7.2). Thus, we have derived the dia- 
grammatic form (7.4) of the fully renormalized two loop scalar QED effective 
Lagrangian in a self-dual background by essentially algebraic manipulations. 
It would be interesting to develop these background field "integration-by- 
parts" rules into a systematic set of rules that might be applied to even 
higher loop order. 

8. Conclusions 

To conclude, I reiterate that from seemingly humble beginnings the 
Heisenberg-Euler result for the one loop effective Lagrangian for spinor QED 
in a constant background field has led to many applications in a wide range 
of branches of modern particle physics and quantum field theory. In some 
sense the Heisenberg-Euler result is analogous to Landau's computation of 
the partition function for nonrelativistic Landau levels which leads to nu- 
merous applications for the magnetic properties of materials. I have only 
skimmed the surface in this review, but hopefully enough groundwork has 
been covered to lead the interested reader directly into the more advanced 
applications. 
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